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Abstract 


Classically, if two highly boosted particles collide head-on, a black hole is expected to form 
whose mass may be inferred from the gravitational radiation emitted during the collision. 
If this occurs at trans-Planckian energies, it should be well described by general relativity. 
Furthermore, if there exist hidden extra dimensions, the fundamental Planck mass may well 
be of the order of the TeV and thus achievable with current or future particle accelerators. 
By modeling the colliding particles as Aichelburg-Sexl shock waves on a flat, D-dimensional 
background, we devise a perturbative framework to compute the space-time metric in the 
future of the collision. Then, a generalisation of Bondi’s formalism is employed to extract 
the gravitational radiation and compute the inelasticity of the collision: the percentage of 
the initial centre-of-mass energy that is radiated away. Using the axial symmetry of the 
problem, we show that this information is encoded in a single function of the transverse 
metric components - the news function. We then unveil a hidden conformal symmetry which 
exists at each order in perturbation theory and thus makes the problem effectively two- 
dimensional. Moreover, it allows for the factorisation of the angular dependence of the 
news function, i.e. the radiation pattern at null infinity, and clarifies the correspondence 
between the perturbative series of the metric and an angular expansion off the collision axis. 
The first-order estimate, or isotropic term, is computed analytically and yields a remarkable 
simple formula for the inelasticity for any D. Higher-order terms, however, require the use 
of a computer for numerical integration. We study the integration domain and compute, 
numerically, the Green’s functions and the sources, thus paving the way for the computation 
of the inelasticity in a future work. 

Keywords: black holes, shock waves, gravitational radiation, extra dimensions, trans- 
Planckian scattering. 
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Resumo 


Classicamente, se duas particulas altamente energeticas colidirem frontalmente, espera-se a 
formagao de um buraco negro cuja massa pode ser inferida a partir da radiagao gravitacional 
emitida durante a colisao. Se isto ocorrer a energias trans-Planckianas, devera ser bem 
descrito pela relatividade geral. Alem disso, se existirem dimensoes extra escondidas, a 
massa de Planck fundamental pode bem ser da ordem do TeV e portanto alcangavel em 
actuals ou futures aceleradores de particulas. Modelando essas particulas como ondas de 
cheque de Aichelburg-Sexl num fundo piano ZD-dimensional, estabelecemos um metodo per¬ 
turbative para calcular a metrica do espago-tempo no future da colisao. Uma generalizagao 
do formalismo de Bondi e entao empregue para extrair a radiagao gravitacional e calcular a 
inelasticidade da colisao: a percentagem da energia inicial do centro-de-massa que e radiada. 
Usando a simetria axial do problema, mostramos que essa informagao esta codificada numa so 
fungao das componentes transversas da metrica - a ‘news function’. De seguida desvendamos 
uma simetria conforme que existe escondida em cada ordem da teoria de perturbagoes e assim 
torna o problema efectivamente bidimensional. Adicionalmente, permite a factorizagao da 
dependencia angular da ‘news function’ e clarifica a correspondencia entre a serie perturbativa 
da metrica e uma expansao angular a partir do eixo. A estimativa de primeira ordem, 
ou o termo isotropico, e calculada analiticamente e produz uma formula simples para a 
inelasticidade em qualquer D. Termos de ordem superior, no entanto, requerem o uso de 
um computador para integragao numerica. Estudamos o dommio de integragao e calculamos, 
numericamente, as fungoes de Green e as fontes, pavimentando assim o caminho para o calculo 
da inelasticidade num trabalho future. 

Palavras-chave: buracos negros, ondas de cheque, radiagao gravitacional, dimensoes extra, 
difusao trans-Planckiana. 
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Chapter 1 


Introduction 


I am going on an adventure! 


Bilbo Baggins 
The Hobbit 


1.1 Why D instead of four? 

The theory of general relativity, presented by Einstein almost a century ago in 1916 [1], 
revolutionised physics by proposing that the gravitational interaction could be described in 
purely geometric terms. It did not take too long for others to realise that this geometrisation 
of apparently non-geometric degrees of freedom could equally describe other fundamental in¬ 
teractions if a higher number of spatial dimensions was invoked: if the world is H-dimensional, 
with D > 4, motion in the ‘extra’ apparently unseen dimensions is not perceived as motion, 
but rather as some other degree of freedom. Since the proposals of Kaluza [2] and Klein 
[3] (preceded by Nordstrom [1]), the converse idea that observed non-motion degrees of 
freedom may be transformed into motion if a higher-dimensional space is considered has 
been a powerful attractor in the quest for a unified description of fundamental interactions. 

In the last few decades, the idea of extra dimensions regained popularity due to its naturalness 
in string theory, which is most naturally formulated in H = 26 (bosonic), D = 10 (perturba¬ 
tive superstring) or H = 11 (non-perturbative superstring or M-theory). Since this framework 
promised to be a fundamental description of all interactions, and in particular a quantum 
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theory of gravity, the study of L>-dimensional gravity was motivated by its appearance as an 
effective low-energy limit of string theory. 

The AdS/CFT correspondence, introduced in 1997 by Maldacena [5], established a duality 
between conformal (quantum) field theories in D dimensions and classical gravity in a (D+l)- 
dimensional anti de Sitter space. This gauge/gravity duality is a realisation of the more 
general holographic principle originally due to t’Hooft [6] and later developed by Susskind 
[7]. In this context, higher-dimensional gravity is motivated by what it can teach us about 
otherwise untreatable problems in gauge theories (both qualitatively and as a computing 
tool). 

However, the study of general relativity regarding the number of space-time dimensions D as a 
parameter can, by itself, be qualitatively and quantitatively informative for the understanding 
of our (apparently) four-dimensional world. Already fifty years ago, Tangherlini [8] showed 
that no stable bound states exist in the D > 4 Schwarzschild solution. A similar argument 
for solutions of the Schrodinger equation provides a reason why atoms or planetary systems 
could not exist in more than four infinite dimensions. This is an interesting lesson, although 
compact dimensions (or brane constructions) can provide a way around this argument. More 
recently, Emparan and Reall [9] demonstrated that the uniqueness theorems for black holes 
do not hold in higher dimensions or, at best, need more assumptions to define uniqueness 
(rod structure, stability, etc). This exemplifies that general relativity has special properties 
in H = 4 which can only be appreciated if H > 4 is considered. 

1.2 The classical trans-Planckian problem 

Relativistic particle collisions are in the realm of quantum field theory and if energies are high 
enough such that gravity becomes relevant they should enter the domain of quantum gravity. 
Moreover, if a black hole forms, as expected in a trans-Planckian head-on collision, non- 
perturbative processes should become relevant and therefore we find ourselves with a hopeless 
problem in non-linear, non-perturbative, strongly time-dependent quantum gravity. However, 
as first argued by’t Hooft [Tn], well above the fundamental Planck scale this process should be 
well described by classical gravity (general relativity), the reason being that the Schwarzschild 
radius for the collision energy becomes much larger than the de Broglie wavelength (for the 
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same energy) or any other interaction scale. 

Thorne’s hoop conjecture HH further tells us that if an amount of energy E is trapped in a 
region of space such that a circular hoop with radius R encloses this matter in all directions, 
a black hole (i.e. an event horizon) is formed if its Schwarzschild radius Rs > R (the classical 
version of this conjecture has been supported by numerical relativity simulations, see for 
example the work by Choptuik and Pretorius on boson stars m)- 

All complex field theoretical interactions will then be cloaked by an event horizon and 
therefore causally disconnected from the exterior. This horizon, in turn, will be sufficiently 
classical if large enough, in the sense that quantum corrections will be small on and outside 
of it. If gravity is the dominant interaction, the total energy of the colliding particles is 
the dominant parameter of the collision. All other constituent details such as gauge charges 
should have a sub-dominant role, i.e. matter does not matter [13]. This idea has been 
verified within numerical relativity in several setups, namely the collision of highly boosted 
black holes, boson stars and self-gravitating fluid spheres naiiaiiei. 

1.2.1 TeV gravity 

The enormity of the Planck mass, mp ~ l/^/GN ~ 10^® GeV, seems to render utopical any 
experimental realisation of this scenario. However, if one admits the possibility of existence 
of extra hidden dimensions, with standard model interactions conhned to a four-dimensional 
brane, the fundamental Planck mass may be well below its effective four-dimensional value. 
Such models were proposed to address what came to be known as the hierarchy problem: the 
relative weakness of gravity by about forty orders of magnitude when compared to the other 
fundamental interactions. Pictorially, it is as if gravity is diluted in the extra dimensions, 
which may be large up to a sub-millimetre scale (models exist with both compact and infinite 
extra spatial dimensions HZi Eg [191120]). 

If the fundamental Planck mass is of the order of the TeV, then high-energy particle colliders 
such as CERN’s Large Hadron Collider (LHC) |2T1[22], or collisions of ultra-high-energy cos¬ 
mic rays with the Earth’s atmosphere [231124] . or even astrophysical black hole environments 
[sgiMlEl], may realise the above scenario with formation and evaporation of microscopic 
black holes. 
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Since 2011, both the CMS [251[2^|3ni and ATLAS [2IIE21E3] collaborations have set bounds 
on such physics beyond the standard model based on the analysis of the 7 TeV LHC data. 
However, these bounds are extremely dependent on regions of parameter space where black 
holes would be produced with masses, at best, close to the unknown Planck scale [Mj- If 
we require the produced objects to be in the semi-classical (and thus computable) regime, 
the cross-sections become negligible at 7 TeV, and only after the upgrade to 13-14 TeV, 
planned to take place in 2015, will the scenario be properly tested. Any improvement in the 
phenomenology of these models is therefore quite timely. 

The two most important quantities for the modeling of black hole production in high-energy 
particle collisions are the critical impact parameter for black hole formation in parton-parton 
scattering and the energy lost in this process, emitted as gravitational radiation. If the latter 
is large and dominates the missing energy, and if it can be calculated with enough precision, 
it could be an important signature for discovery or exclusion. The event generators used at 
the LHC to look for signatures of black hole production and evaporation, charybdis 2 |35j 
and BLACKMAX [36], are very sensitive to these two parameters. 


1.3 Current estimates of the inelasticity 

In the highly trans-Planckian limit, the colliding particles are greatly boosted, traveling very 
close to the speed of light. This has motivated the study of gravitational shock wave collisions 
as a model for the gravitational fields of highly boosted particles. 

1.3.1 Apparent horizon bounds 

The gravitational field of an ultra-relativistic particle of energy E is obtainable from a boost of 
the Schwarzschild metric. As the boost increases, the gravitational field becomes increasingly 
Lorentz contracted and in the limit in which the velocity goes to c (but keeping the energy 
fixed) the gravitational field (i.e. tidal forces, described by the Riemann tensor) becomes 
planar and has support only on a null surface. This shock wave is described by the Aichelburg- 
Sexl metric m- Due to their flatness outside this null surface, it is possible to superimpose 
two oppositely moving waves and the geometry, as an exact solution of general relativity, is 
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completely known everywhere except in the future light cone of the collision. 

Strikingly, such knowledge is enough to actually show the existence of an apparent horizon 
for this geometry and thus provide strong evidence that a black hole (i.e. an event horizon) 
forms in the collision. Thus, if cosmic censorship holds, and since the sections of the event 
horizon must lie outside the apparent horizon, the size of the latter yields a lower bound on 
the size of the black hole. An energy balance argument then provides an upper bound on 
the inelasticity e, i.e. the percentage of initial centre-of-mass energy that is radiated away as 
gravitational radiation. 


Penrose pioneered this computation in Z) = 4, and by Ending an apparent horizon on the 
past light cone of the collision concluded that no more than 1 — l/-v/2, or about 29.3%, of 
the energy was lost into gravitational waves. His method was later generalised to arbitrary 
D dimensions by Eardley and Giddings [38], who obtained 


EAH < 1 


1 7D — 2 

2 V 2 Hd-J 


(1.3.1) 


where 0^ is the volume of the unit n-sphere. Note that this bound increases monotonically 
approaching 50% as H —?• oo. Later, Yoshino and Rychkov found an apparent horizon in 
the future light cone |39|. Their analysis coincided with the bound in Eq. (1.3.1) for head- 


on collisions, but the critical impact parameters they obtained were larger than previous 
estimates 


1.3.2 A perturbative approach 

Instead of computing a bound one may decide to compute the precise inelasticity by solving 
Einstein’s equations in the future of the collision. This is a tour de force. A method, first 
developed by D’Eath and Chapman and later by D’Eath and Payne HHIEiQiii , is to set 
up a perturbative approach: considering the collision in a highly boosted frame, one shock 
becomes much stronger than the other and the latter can be considered as a perturbation 
of the former. This is how they justified, conceptually, the perturbative expansion. In four 
dimensions, they obtained a value of 25% for the inelasticity at first order in perturbation 
theory, in agreement with Smarr’s Zero Frequency Limit [Hj. This was originally thought 
to be the exact value, but Payne |35| showed that Smarr’s formula is in fact a linearised 
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approximation valid only when the gravitational radiation is weak, and cannot predict the 
strong-field radiation generated by fully non-linear gravitational interactions. Second-order 
perturbation theory lowered the inelasticity to 16.4%. This is within the range 14 ± 3% 
obtained by considering the collision of highly boosted black holes in numerical relativity 
m, which confers further validity to the approach. Moreover, these values are smaller than 
the apparent horizon bound. If less energy is lost into gravitational radiation, the final black 
hole is more massive and hence more consistent with the semi-classical analysis used for the 
potentially observable Hawking radiation. 


1.4 Collision of gravitational shock waves in D dimensions 

Phenomenologically interesting TeV gravity models occur in dimensions greater than six. 
Therefore, it would be helpful to extend the above-mentioned methods to higher dimensions. 
Numerical relativity for higher-dimensional space-times has seen an increasing amount of 
development in recent years uni HZ], but so far black hole collisions have only been successfully 
computed at low energies [381 US] • 

Extending the method of D’Eath and Payne to higher generic H is a demanding task involving 
analytical and numerical methods. This was the goal set out by the Gr@v group at the 
University of Aveiro in 2011 just prior to the author joining as a PhD student. The remainder 
of this thesis will describe our findings and conclusions so far. Its structure is as follows. 

We start in Chapter by studying the geometrical properties of Aichelburg-Sexl shock waves 
and setting up a space-time where two such waves collide head-on. 

Then, in Chapter we devise a perturbative strategy to compute the metric in the future of 
the collision; using the axial symmetry of the problem, we reduce it to three dimensions and 
write the general solution for the metric at each order in perturbation theory. 

Chapter is about gravitational radiation and how to extract it in our problem, using a 
generalisation of the Bondi formalism. 

Then, in Chapter we unveil a hidden symmetry of the problem which allows for the 
factorisation of the angular dependence of the radiation pattern at null infinity. 


1.4. COLLISION OF GRAVITATIONAL SHOCK WAVES IN D DIMENSIONS 
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With this, we compute analytically the first-order estimate for the inelasticity in Chapter 

Next, in Chapter we use this hidden symmetry to reduce the problem to two-dimensions 
and prepare it for numerical evaluation in a computer; we find characteristic coordinates and 
produce a conformal diagram of the effective two-dimensional space-time. 

In Chapter we pave the way for the second-order calculation by computing numerically the 
Green’s functions and the sources, and discussing the asymptotic behaviour of these functions 
at the boundaries of the integration domain. 

We conclude in Chapter with some final remarks and an outlook on the second-order 
calculation, which we hope to finish in a future work. 

Many technical details are given in Appendices. 



Chapter 2 


Kinematics of shock wave collisions 


If you can’t explain it simply, you don’t 
understand it well enough. 


Albert Einstein 


2.1 The Aichelburg-Sexl shock wave 


The Aichelburg-Sexl (AS) metric is an exact solution of general relativity representing the 
gravitational field of a point-like particle moving at the speed of light. It was first obtained 
in ID = 4 by Aichelburg and Sexl |37j by boosting the Schwarzschild solution and taking the 
speed of light limit while keeping the energy finite. 


For higher dimensions, the starting point is the ID-dimensional Schwarzschild-Tangherlini 
solution, i.e. a spherically-symmetric, static black hole with mass M in ID dimensions, 

-1 


ds^ = - l- 


IGttGdM 


1 


(ID - 2)Q.d- 2 r 


D-3 


dr+1- 


IGttGdM 


1 


(ID — 2)0£)_2 r 


D-3 


dr^ + r^dQjj_2 , 


( 2 . 1 . 1 ) 

where dQ‘jy _2 and ^d -2 are the line element and volume of the unit (ID — 2)-sphere and Go 
is the ID-dimensional Newton constant. Before taking the boost, it is convenient to rewrite 
this in isotropic coordinates by changing to a new radial coordinate R, 


r = R{1 + Z)o-3 , 





WttGdM 
(ID — 2)fD25_2 


( 2 . 1 . 2 ) 
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Here r* is the Schwarzschild radius of the black hole and Z is simply a short-hand. Then Eq. 


( |2.1.1[ ) becomes ( R = \lY^f=i 


ds^ = — 


1 - Z 

TTz 


D-l 


dt^ + {1 + Z) d-3 ^ dx'f . 


(2.1.3) 


i=l 


Next, we choose a spatial direction (say along ^ = z) and perform a boost with velocity 
i.e. 

t —>■ 7 (t — vz), z —>■ 7 ( 2 ; — vt ), (2.1.4) 

where 7 = (1 — u ^)~2 is the Lorentz factor. 

Meanwhile the transverse coordinates xt = (xi,... ,xd- 2 ) are unaffected by the boost. The 
quantity Z now reads 


Z = 


A-kGdE 


1 


{D 2)H£)_2 ^ 2 


(2.1.5) 


where E = 7 M is the energy and p = radius in the transverse plane. 

In the limit u —1, Z —)• 0 and the metric can be expanded as 


ds^ = —dP + dz^ + dxr + 4 


.D-2 


D-3 


'y‘^Z{dt — vdz)"^ + ... , 


( 2 . 1 . 6 ) 


where ... denotes terms that vanish in the limit. 


Fig. 2.1 provides an illustration of this procedure: in a frame where the black hole is moving 
with velocity v, the gravitational field is Lorentz contracted along 2 : and the curvature becomes 
increasingly concentrated on a plane perpendicular to the motion; transverse directions are 


not affected. Indeed, in the limit u —)• 1 the term proportional to Z in Eq. (2.1.6) goes to 


zero off the plane z = t (which is moving at the speed of light) and diverges on it. 

In Appendix |V we show that the end result is simply flat Minkowski space-time plus a Dirac 
delta function with a radial profile on the transverse moving plane, 


= —dudv -|- d/ 9 ^ -b p‘^dLljj_^ -b K^{p)S{u)du "^, 


(2.1.7) 


where k = 8 ttGdE/Qd -3 and u = t — z,v = t + z are null coordinates. The profile function 
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Figure 2.1: Schematic representation of the qualitative behaviour of the gravitational field 
(Kretschmann scalar —\R^'^°‘^Rfj.^a/ 3 \) of a black hole subject to a boost (horizontal 
directions represent the isotropic spatial coordinates Xi) at some fixed time. As the 
boost increases, the field becomes increasingly concentrated on the transverse plane 
containing the particle. From m- 


$ depends only on p and takes the forn:[^ 

$(p) = < 


-21n(p) , 
2 


D = A 
D>A 


( 2 . 1 . 8 ) 


{D - A)pD-^ ’ 

Clearly, a shock wave moving in the opposite direction with the same energy is obtained by 
replacing 2 ; •(-)> — 2 ; or, equivalently, exchanging u and v. 


2.2 Properties of the Aichelburg-Sexl solution 


The Aichelburg-Sexl metric, Eq. (2.1.7), has the following properties: 


• axial symmetry: it is invariant under rotations {(NId-s) on the transverse plane; 

• n-translational symmetry: the metric components do not depend on v, 

• transformation under boosts: a further boost of the solution moving along -I- 2 ; with 
velocity v = tanh /3 amounts to a rescaling of the energy parameter E —>■ e^E, as 
expected from the transformation law of the H-momentum of a null particle. 


It is also instructive to look at the Riemann tensor in order to better understand the 


gravitational held. One can either compute it directly from Eq. (2.1.7) or from Eq. (2.1.1) 


*In D = 4 the profile contains an arbitrary length scale, ^{p) = —2 ln(p/t), but since all physical quantities 


will be independent from it, we set (. = 1. 
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Figure 2.2: Profile function 4>(/o) of the gravitational field on the transverse plane, for = 5 (left) 
and D = A (right). The plane is shown in perspective and in the inset smaller figures 
we plot $ as a function of a transverse direction Xi. From isl¬ 


and then applying the boost and taking the limit. The only independent non-vanishing 
component is 




UlUJ 






( 2 . 2 . 1 ) 


where is the Laplacian on the transverse plane. In Eq. (2.2.1) it is manifest that the 
gravitational field only has support on the plane u = 0, decaying away from its centre 
as an inverse power of p. This can be better visualised in Fig. |2.2| showing plots of the 
profile function on the transverse plane. For D = 5 {D > A hi general), 4>(/3) decays for 
large p, whereas in H = 4 it goes to —oo. This special case does not contradict the above 
statements on the qualitative behaviour of the gravitational field, since $ is not a gauge- 
invariant quantity. Any gauge-invariant object contains at least one derivative like <!>' in 


Eq. (2.2.1), so the gravitational field does indeed decay for large p. 


From the Riemann tensor one constructs the Einstein tensor or, equivalently, the energy- 
momentum tensor, whose only non-vanishing component reads 


T — 

im. — 


E 


6{u)V^<I> = ES{u)5^^-^\xi). 


( 2 . 2 . 2 ) 
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It is now clear that this solution describes the gravitational field of a point-like particle moving 
at the speed of light with energy E. Furthermore, the impulsive nature of this shock wave, 
encoded in the 5{u) and distributions, can be seen as an idealisation of a more 

realistic distribution replacing the delta function by a smeared out form, as would be the case 
if the matter distribution had some spatial extension. As long as this source is sufficiently 
localised, we expect the gravitational held away from the centre to be insensitive to these 
details. 


2.3 Geometric optics 


We will now look at the null geodesics of this space-time in order to understand the grav¬ 
itational dehection and redshift suffered by null test particles. Later, when considering the 
collision of two AS shock waves, such null rays will also be useful to dehne the causal structure 
of the post-collision space-time. 

The parametric equations of null geodesics can be given as 


u 

V 


u, 

Vo + K^{po)9{u) + u 

i 

Xq + x'qU + ^Vi^(po)u9(u), 


(2.3.1) 


where the coordinate u is itself an affine parameter, vq, Xq and x'q are integration constants 
and 9(u) is the Heaviside step function. In these coordinates (known as Brinkmann coor¬ 
dinates m) geodesics and their tangent vectors are discontinuous across the shock wave. 
This problem is solved by identifying the parameters of null geodesics which are incident 
perpendicularly to the shock plane (x'q = 0) with new coordinates (h, i;,x*). 


u 

V 


u, 

V -|- k9{u) 



X* + '^Vi^(x)u9(u) 


P 


p(l+'^^'^ 

4*a 


(2.3.2) 


where <I> and are evaluated at p and a = 1... D — 3, are the angles on the (T>—3)-sphere. 
These are known as Rosen coordinates |38j |52j. 
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From inspecting Eqs. (2.3.2) we can draw the following conclusions: 


1. Redshift: null rays incident perpendicularly to the plane of the shock wave {u = u = 0) 
suffer a discontinuous jump in v, which is always positive for D > 4 but becomes 
negative for large p in D = 4. 

2. Focusing: they are further deflected by an angle obeying tana = and focus, 

generating a caustic, at 


p = 0, 


V = V + K < 


1 — In(KM) 


D-2 


D-4 


(ku) 


D-4 

D-2 


D = 4 
D>4 


(2.3.3) 


This focusing increases (decreases) with D for short (long) distances, as expected from 
the behaviour of the gravitational force. 


These properties are illustrated in Fig. 2.3, where the small red arrows represent the scattered 
trajectories of a plane of null rays and the big blue arrow represents an AS shock wave moving 
to the right. The rays propagate freely (i.e. in a straight line) before crossing the shock (left 
panel). They are scattered instantaneously at t = 0 when the two planes meet and thereafter 
continue on another straight line. 


The two effects described above are clearly seen in the middle and right panels. In the middle 
one, the tangent vectors are bent inwards towards the axis, as expected from the gravitational 
attraction of the source at the centre - this is the focusing effect. Moreover, the farther the 
rays are from the axis, the lesser they are bent. Later (right panel), we see there is an 
increasingly circular envelope of rays around the scattering centre, and another outermost 
one which asymptotes the initial wavefront (dashed line). 

The redshift effect can be inferred from the middle panel: we see that there are rays which 
have not yet emerged from the central region of the shock wave - they get stuck for some 
time until emerging later in their bent trajectory. Again, this effect is stronger the closer the 
rays are to the centre of the shock wave. 
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Figure 2.3: Representation of a test wavefront of null rays (small red arrows) scattering through the 
transverse plane of the AS shock wave travelling along the positive z direction (big blue 
arrow), for D = 5: i) before scattering (left), ii) a small time after scattering (middle) 
and iii) a later time after scattering (right). The dashed red line (right) represents the 
plane wavefront the rays would have followed if no scattering had occurred. From m- 


2.4 Superposition and the apparent horizon 


In the previous two sections we showed that an AS shock wave with support on u = 0 
represents a massless particle moving along z = t at the speed of light. Moreover, we saw 
that its gravitational field is restricted to the transverse plane. In particular, an incident ray 
is not affected before it crosses the shock. This is simply due to causality: since no signal 
can travel faster than the speed of light, the incident ray has no means of “knowing” what is 
coming before it hits the shock and thus experiences no gravitational held. 

Therefore, if we superpose another shock wave travelling in the opposite direction, i.e. along 
z = —t, the space-time metric is simply the superposition of the two metrics everywhere but 
in the future of the collision. 


In Rosen coordinates, the line element (2.1.7) becomes 


ds^ = —dudv -b ( 1 + 


nu9{u) 




dp^ + 


1 ^ug(u) ^, 

2p 




(2.4.1) 


The geometry for an identical shock wave traveling in the opposite direction is obtained by 
exchanging u ^ v. If the shocks have different energy parameters v and A, the superposed 
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Figure 2.4: Space-time diagram of the collision (left) and apparent horizon upper bound on the 
inelasticity (right) as a function of D. From 


metric reads 


ds'^ = —dudv + 


^ vuOiu) 
1 H- 


+ 1 + 


>^ve{vi „ 


_ 1 


dp^ 


+f 


1 

2-p 


, ^ \v6iv) , 

+ ( 1 + 1 - 1 


dn"^ 


D-3 ) 


(2.4.2) 


and is valid everywhere except in the future light cone of u = u = 0. 


Fig. 2.4 is a space-time diagram, in Rosen coordinates, where the transverse dimensions have 
been suppressed. The two shock waves are represented by the u, v axes where they have 
support. Regions I, II and III are flat. Region I is the region in-between the shocks, before 
the collision, whereas II and III are, respectively, the regions behind the u and v shocks. 

The apparent horizon mentioned in Chapter was found, first by Penrose in D = 4 and later 
by Eardley and Giddings [38] in higher D, on the union of the two null surfaces ft = OAu < 0 


and u = 0 A u < 0. Their upper bound on the inelasticity, Eq. (1.3.1), is plotted in the right 


panel of Fig. 2.4 where we see that it increases monotonically with D, approaching 50% in 
the limit D —)• oo. 


The existence of an apparent horizon, which traps the colliding particles, further supports the 
independence of the results on the details of the localised sources. In particular, their point¬ 
like nature (as opposed to a small but finite volume) is not expected to affect the formation 
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of a black hole in the trans-Planckian limit. 

In the next chapter we shall devise a perturbative framework to compute the metric in the 
future of the collision, i.e. in region IV. 




Chapter 3 


Dynamics: a perturbative approach 


Physicists like to think that all you have to 
do is say, ‘These are the conditions, now 
what happens next?’ 


Richard Feynman 
The Character of Physical Law 


3.1 The collision in a boosted frame 


The original argument of D’Eath and Payne consisted in looking at the collision in a boosted 
frame, in which one of the shocks appears much stronger than the other. This allowed them 
to consider the latter as a perturbation of the former, and set up an iterative process of 
solving Einstein’s equations order by order in perturbation theory. 


Starting from the Rosen form of the metric of one shock wave, Eq. (2.4.1), a boost with 
velocity v in the fuz direction amounts to a scaling of the energy parameter. 




/3 = tanh v . 


(3.1.1) 


If two identical but oppositely travelling shocks are superposed, this change of reference frame 
makes them appear to have different energy parameters. The metric, in this case, would be 


given (except in the future of the collision) by Eq. (2.4.2), with 

h' = e^K , A = e~^K . 


(3.1.2) 
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In this way, the i^-shock (which we call the strong shock) carries much more energy than the 
A-shock (which we call the weak shock), since 

- = e^P = -^oo. (3.1.3) 

A 1 — V 

Thus one may face the shock wave travelling in the — ^ direction as a small perturbation of 
the shock wave travelling in the +z direction. Since the geometry of the latter is flat for u > 0 
(regions II and IV), we can make a perturbative expansion of Einstein’s equations around 
flat space, in order to compute the metric in the future of the collision (h, u > 0, i.e. region 
IV). 


Fig. 3.1 illustrates the effect of the boost in the pattern of gravitational radiation. In the 
boosted frame, the whole z < 0 hemisphere maps to a region around 0 = vr in the centre- 
of-mass frame. Thus that is the region where perturbation theory is expected to be valid, 
with higher orders needed to extrapolate off the axis. This statement will be made clear 


in Section 5.1.3 where we shall derive an exact correspondence between a perturbative 
expansion of the metric (in the sense here described), and an angular expansion of the 
radiation pattern at null infinity. 


3.2 The collision in the centre-of-mass frame 


The independence of physics from the reference frame or the coordinates used is one of 
the pillars of general relativity. Therefore, one must be able to construct the perturbative 
approach of the previous section in the centre-of-mass frame in which the two colliding shock 
waves have the same energy parameter. 


We start by rewriting the superposed metric, Eq. (2.4.1), for identical shocks {u = X = k), 


in a form which resembles a perturbation of flat space. 


ds^ = —dudv + Sijdx^dx^ + 


( K^' 


/\ 2 


+ - 


2 {ue{u) + ve{v)) Aij + 

[if6{u) -I- v‘^6{v)) [{D - 3)6ij - {D - 4)Ay) 


dx^dx^ , (3.2.1) 


V 2p 

where Ajj = Sij — [D — 2)rjrj is a traceless tensor on the transverse plane and Tj = Xi/p are 


angular factors. The first two terms in Eq. (3.2.1) are simply the Minkowski line element 
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Figure 3.1: Space-time diagram of the collision in the centre-of-mass (left panel) and boosted (right 
panel) frames. The transverse a;* directions are orthogonal to the t and z axes. In 
the centre-of-mass frame, the radiation is symmetric under z —>■ —z. However, only the 
radiation that propagates along a small solid angle around 0 = tt in that frame (red area, 
delimited by the polar angle tand = — 1) is still propagating in the —z direction 

in the boosted frame. All the radiation in this solid angle is redshifted in that frame. 
From JS^y . 


and those inside the square brackets constitute a small perturbation in regions where 

P»{ku)t^ a p»{k,v)t^. (3.2.2) 

These two conditions highlight the main shortcoming of Rosen coordinates: in regions II and 
III the space-time is exactly flat, but the metric is not in standard Minkowski form, and 
actually depends on the energy k of the shock waves. This can be partially fixed by returning 


to Brinkmann coordinates through Eq. (2.3.2). Since this transformation is adapted to the 


shock wave with support on n = u = 0, the z ■H' —z symmetry of the problem is no longer 
apparent. 


By splitting the u-shock from the n-shock, Eq. (3.2.1) can be written as 


ds'^ = —dudv -b dijdx^dx^ -)- K^{p)6{u)du^+ 

{—2h{u, V, p)Aij + h{u, V, p)^ [{D — 3)6ij — {D — 4)Aij) | dxPdx^ , (3.2.3) 

where the barred coordinates in the second line are to be understood as functions of the 


unbarred coordinates through Eq. (2.3.2), and 
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Figure 3.2: Right: 3D space-time diagram showing the {t,z,p = x) axes, the numbered regions (in 
red), the generators of the u-shock (in blue) and the collision surface (in green). From 

Left: Top diagram shows the surface defined by the null generators of the u-shock as 
they scatter through the u-shock. Bottom one shows the causal surfaces (u = 0 A u > 0 
and u = 0 A It > 0) defining the future light cone of the collision (region IV). From 15d| /. 


As long as h{u,v,p) « 1, the tt-shock now looks like a small perturbation on a u-shock 
background (which is itself flat everywhere except on u = 0). This is analogous to the picture 


of Section 3.1 where we saw that, in a boosted frame, one shock looks much stronger than 
the other. 


At this point, it is instructive to revisit Section 2.3 and repeat the analysis of the null rays 
in the new adapted Brinkmann coordinates. First note that, from the n-shock point of view, 
the generators of the u-shock (i.e. the null rays travelling on u = 0) will scatter through the 
u-shock exactly as the incident ‘test’ rays of Fig. |2.3[ Thus the outermost envelope of rays. 


better seen in the right panel of Fig. 2.3, defines the causal boundary between the curved 
region (IV) and the flat region (II). 
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The right panel of Fig. 3.2 represents a space-time diagram of the collision in u-adapted 
Brinkmann coordinates. The generators of the tt-shock travel along u = 0 without any 
discontinuity. The generators of the u-shock, however, travel along u = 0 while u < 0 (blue 
arrows) but suffer a discontinuous jump to the green line upon hitting the ‘strong’ shock at 
u = 0. The green line is the collision surface defined by ii = h = 0, or 


u = 0 A v = <L(p}. (3.2.5) 

After that they focus towards the axis and cross, forming a caustic (red line). Region I’ 
represents a patch of flat space below region IV in u > 0 where the null rays have not yet 
crossed and thus space-time is still flat. 

The 3D diagrams on the left complement the visualisation. The top one shows the surface 
formed by the generators of the u-shock, together with the world-line of a time-like observer 
for comparison. The bottom one represents the boundary of the curved region (IV). 


3.3 An initial value problem 


In this section, we shall use the results of the previous one to set up an initial value problem 


for the dynamical equations that determine the metric in region IV. From Eq. (3.2.1) we see 


that onu = 0Au>0 (blue surface on the bottom left diagram of Fig. 3.2), the metric has a 
standard Minkowski form, 

g^iy(u > 0,v = 0,x') = . (3.3.1) 

On the other hand, on u = O'*' (red surface). 


g/,uiu = 0+, u, X*) = rif,u + + /ijfj , (3.3.2) 

where the k in is defined by the power of h{u = 0^,u,x*) it contains. In natural units 
where k = iQ 

Ai’ = -{0-3)4>'ft(t>.P)U. = (3.3.3) 

= -2h{v, p)Aij , hf^ = h{v, pf {{D - 3)6ij - {D - 4)Aij) , 

* K has dimensions of [Lenght]^“'^. 









46 


CHAPTER 3. DYNAMICS: A PERTURBATIVE APPROACH 


where 


h(v, p) = h{0, v,p) = ~{v - <I>{p))9{v - <^{p)) 

2p 


(3.3.4) 


These initial conditions, Eqs. (3.3.1) and (3.3.2), have the form of a perturbation of a flat 


Minkowski background which is exact at second order (on u = 0^). 

By adopting the reference frame of the u-shock, we have successfully encoded the information 


on the scattering of the u-shock in Eq. (3.3.2). It is important to stress that no approximation 


has yet been made - all we did was chose a convenient frame for the superposition. 

To find the metric in the future of the collision (region IV) we must solve Einstein’s equations 


in vacuum, Rfj,^ = 0, subject to initial conditions on u = 0^. The form of Eq. (3.3.2) suggests 


the use of perturbation theory, in which one obtains each new term of an infinite series in an 
iterative process. However, and are only small in regions where 


\h{v,p)\ <C1 ^ g{p) = 




2p 


4>'(p) 


> u. 


(3.3.5) 


This function g{p) is plotted in Eig. 3.3, where we see it grows very fast for large p. Combined 


with the fact that rays incident close to the axis are trapped inside the apparent horizon 
(which is at p = yO = 1 ), this suggests that the outside region of the resulting black hole 
might be well described by perturbation theory. 


Eurthermore, we saw in Section 3.1 that perturbation theory should be valid in the region 
close to 0 = TT, i.e. close to the —z axis. Indeed, due to focusing, the rays reaching such 


an observer do come from the far field region, as shown in the right panel of Eig. 3.3 and 
in Eig. |3.4[ This is the true meaning of perturbation theory in this problem; it gives an 
angular expansion of the metric in the future of the collision, with higher orders needed to 
extrapolate further off the axis or 0 = tt. However, a quantitative restatement of this idea 
will have to wait until Chapter 


3.3.1 The dynamical equations 

The standard way of applying perturbation theory in general relativity starts with an ansatz 
for the metric, 

OO 

9^lU = ^ ^ ) 

k=l 


(3.3.6) 
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9{\xi\) 



Xi 



Figure 3.3: Left: Function g{p) which determines the validity of the perturbative approximation. 

Middle: Diagrammatic representation of the u-shock null generators for large p in D = 5. 
Right: Spatial slice (in D = 5) showing the observation region behind the collision event 
reached by such rays with large p, where perturbation theory is expected to yield a good 
result. From m- 


(k) 

where each term is assumed to be sufficiently smaller than the previous ones such that the 
series converges or a consistent truncation can be made. Then Einstein’s equations produce 
a tower of linear tensorial equations, one for each order, where the source for each h)iP is 
given by an effective energy-momentum tensor built out of lower order terms. 


Each of these equations generally couples different components, but it is well known that 
coordinates can be chosen such that they decouple. Defining the trace-reversed metric 
perturbation 

, (3.3.7) 

the de Bonder or Lorentz gauge is obtained by imposing D conditions 


= 0 . 


(3.3.8) 


Then the field equations become a set of wave equations for each component. 




‘‘ap 


(3.3.9) 


where □ = is the D-dimensional d’Alembertian operator and the source on the 

right-hand side can be computed explicitly order by order. 


Here and throughout this thesis, the superscript index {k) denotes the order of the space- 
time perturbation theory used (i.e. needed) to compute the given object. In particular 
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Figure 3.4: Illustration of the spatial trajectories of the weak shock null generators, exhibiting their 
focusing after u = 0. A point on the axis {p = 0) will be hit, simultaneously, by an infinite 
number of null generators (two of each are in the represented section). This occurs at 
the caustic and marks the beginning of the curved region of space-time. Therefore, 
the continuation of the rays beyond the axis, represented by a dotted line, is merely 
illustrative. For points outside the axis, such as point P, this diagram suggests an initial 
radiation signal, associated with ray 1, followed by a later burst, associated with ray 
2. The closer the observation point V is to the axis, the farther those rays come from. 
From 


observe that, for objects which are not linear in the metric, this does not correspond to the 


perturbation theory order of the object itself. For example, in Eq. 3.3.9 at 0{k), the right- 

(k—l) 

hand side is denoted ^ because, being quadratic in the metric, it is built out of lower 
order metric perturbations. 


3.4 Formal solution 


The formal solution to Eq. (3.3.9) can be written using the Green’s function for the □ 
operator. In Appendix we define it as the solution to the equation 


□G(x'^) = <5(x^), 


(3.4.1) 
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and show that, in our coordinates = {u, v, x*}, it is given bsQ 


-^ Hx), X =. 

dvr 2 


(3.4.2) 


Then, using Theorem 6.3.1 of m, the metric perturbations in de Bonder gauge are given by 


hl^J{x) = F.P. I d^x'G{x-x') \46{u')d,>h^^X^') + Tlt-^\^') 


"'{lU ' 


' flU 


(3.4.3) 


where F.P. denotes the hnite part of the integral and the source only has support in u > 0. 


Using Eq. (3.4.1), it is straightforward to check that Eq. (3.4.3) obeys Eqs. (3.3.2) and (3.3.9) 


(see Appendix (B.1.1)). 


The first term inside the square brackets of Eq. (3.4.3), which we shall call the surface 


term, simply propagates the initial data on u = O"*" and makes the solution compliant with 


Eq. (3.3.2). The second one, which we shall call the volume term, encodes the non-linearities 


resulting from the interference with the background radiation sources generated by lower 
order terms. 

Note that at first order there is no source: the solution is simply the propagation of the initial 
conditions on a flat background. Moreover, since the initial data are exact at second order, 
for k > 3 there are only volume terms. This might be the reason why the second order result 
of D’Eath and Payne in B = 4 is already in good agreement with the predictions of numerical 
relativity (see Chapter Q. 


Fig. 3.5 is a space-time diagram showing the past light cone of an observer P after the 
collision. This is defined by x = 0 and its intersection with the u = 0 plane is a parabola. 


lx* — x'*P 


V = V — 


(3.4.4) 


The initial data on tt = 0 have support above the green line of Fig. 3.2 (where the collision 
takes place). 


denotes the n-th derivative of the delta function. See Appendix B for the meaning of non-integer n. 
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Figure 3.5: The dashed blue parabola is the intersection of the past light cone of a space-time event 
V{u, V, x), to the future of u = 0 = ?;, with the null hypersurface u = 0. From M- 

3.4.1 Gauge fixing 


The initial data, Eqs. (3.3.3), are not in de Bonder gauge. Therefore, before inserting them 


in the formal solution, Eq. (3.4.3), we need to make a coordinate transformation 

OO 


(3.4.5) 


k=l 


where will be determined order by order through the de Bonder gauge condition, Eq. (3.3.8). 
In what follows we derive a condition to be imposed on the initial data, such that the de 


Bonder gauge is preserved by the evolution equations. From Eq. (3.3.9), 




(3.4.6) 


where, as before, the bar denotes trace-reversed quantities. The conservation of the effective 
energy-momentum tensor was checked explicitly up to third order, but it must also hold at 
higher orders since it is a consequence of the symmetries of the (flat vacuum) background. 
The formal solution to the de Bonder gauge condition is then 

hg)-‘^(x) = 4 y G{x - x')d(u')5^'hjfj’^(x'), (3.4.7) 

where the hnite part prescription was left implicit. Thus we conclude that 



^= 0 + 


= 0 , 


(3.4.8) 


is a sufficient condition for the de Bonder gauge to be maintained for all u. 
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3.5 Reduction to three dimensions 


The general solution of Eq. (3.4.3) does not take into account the axial symmetry of the 


problem (rotations around the z axis). Indeed, if we work with spherical coordinates on the 
(D — 2)-dimensional transverse space, we should be able to integrate out the D — 3 angles 
(/>“, and be left with an integral in p = x^Xi only. 

In the coordinates we have been using so far, = {u, v, x*}, the metric has the generic form 


ds'^ = ds'jriat + huudu^ + 2huidx^du + hwdv^ + 2huvdudv + 2hmdx^dv + hijdx^dx^ . (3.5.1) 


Dehning a basis of tensors in that plane, 

5ij, Ti = ^, Aij = 6ij - {D - 2)TiTj, (3.5.2) 

where we have chosen the last one to be traceless, we can decompose the metric perturbations 
hij into seven functions of (n, v, p), here denoted A, B, C, E, F, G, H, in the following way: 

Ku^A = + ^(2) + ... Ki^BT, = (S(1) + 5(2) + .. .)ri 

Kv^C = C7(1) + C(2) + ... Ki = FTi = (F(i) + F(2) + .. .)ri (3.5.3) 

hy^ = G = G(i) + G(2) + ... h^J = E Aij + H 6ij 

= (5(1) + .. .)\, + (5(1) + .. .)dij . 


(k—1) 

A similar decomposition can be applied to the sources T^i, ^. For example, for the transverse 
components, 

Rik-i) ^ r(f-')(u,u,p)(I.,+r(f-')(u,u,p)A,,. (3.5.4) 

Now let F denote generically any of those functions and let S be its associated source, 
including both volume and surface terms, 

S^’‘\u,v, p) = T^^~^\u,v, p) + 4:S{u)dvF^'^\o,v, p) . (3.5.5) 


In Section B.2 of Appendix [B] we show that Eq. (3.4.3) implies 

F{u,v,p)= / du' / dv' / dp'Gm{u — u',v — v'■, p, p')S{u p'), 


(3.5.6) 
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where the new Green’s function Gm depends on the rank m = {0,1,2} of the invariant tensor 
multiplying F in the respective component of h^,^. For example, m = 0 for A and H, m = 1 


for B and m = 2 for E as can be observed in Eqs. (3.5.3). 

It is given by 

D-4 

Gm{u-u',v-v';p,p') = 


where 


X* = 


p^ + p'^ — {u — u'){v — P) 


and the one-dimensional functions /m’” are defined as 

r-l 




, ^ D-2 

(2vr)^ 


(3.5.7) 

(3.5.8) 

(3.5.9) 


with Am = {1 ,-{P - 3), {D - 1)(F) - 3)}. 


Thus, the problem of finding the post-collision metric is now formally solved. The extraction 
of physically relevant information from that solution will be the subject of the next chapter. 





Chapter 4 


Extraction of gravitational 
radiation 


energy • the property of matter and ra¬ 
diation which is manifest as a capacity to 
perform work; from the greek energeia. 


Oxford English Dictionary 


In Chapter we analysed the geometric properties of shock waves and in Chapter we 
developed a perturbative framework which allowed us to obtain a formal expression for the 
metric in the future of the collision. It is now time for us to see how this can be used 
to extract physically meaningful information, namely the energy that is radiated away by 
gravitational waves. If the end product of the collision is a black hole at rest, its mass must 
be the difference between the initial energy of the colliding particles and the total energy 
radiated. We shall see that the power or flux of gravitational energy that reaches a sphere at 
‘infinity’ is characterised by a single function of a retarded time r and angular coordinate 9, 
called the news function, which is then used to derive a formula for the inelasticity. 


We start in Section 4.1 by discussing the concept of energy in asymptotically flat space-times 


and introducing the Bondi mass, which is then detailed in Section 4.2 for arbitrary D di¬ 
mensions. Section [4.3| establishes the relationship between de Bonder and Bondi coordinates. 


Due to its technical character, a summary is provided in the end, Section 4.3.1, to which a 
reader may safely jump without loss of comprehensibility. 
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4.1 Mass and energy in asymptotically flat space-times 


The notion of mass, or energy, in general relativity has always been a delicate topic due to the 
difficulty in defining it precisely in generic terms. In particular, there is no local definition of 
the energy density of the gravitational field (as there is, for instance, for the electromagnetic 
field). However, satisfactory dehnitions of total energy and radiated energy have been found 
for isolated systems, for which the behaviour of the gravitational field at large distances away 
from the source can provide a notion of ‘total gravitational mass’. 

Ideally isolated systems live in what are known as asymptotically flat space-times. Several 
definitions exist, with varying degrees of mathematical sophistication (see for example P]), 
but they all rely on some notion of ‘infinity’ and on appropriate fall-off behaviour of the 
metric components (in suitable coordinates). 

For stationary, asymptotically flat space-times (which are vacuum near infinity), the time 
translation symmetry implies there is an associated conserved quantity, called the Komar 
mass P], which is given by an integral over a distant sphere of the corresponding Killing 
vector field. In the non-stationary case, alternative definitions must be found. There are two 
main options. 

The Arnowitt-Deser-Misner (ADM) energy-momentum was motivated by the Hamiltonian 
formulation of general relativity [58], and gives a well-defined notion of energy/mass and 
momentum at spatial infinity. It is usually interpreted as the total energy available in the 
space-time. In particular, it cannot differentiate between rest mass and energy stored in 
gravitational waves, which makes it inappropriate for our purposes. 


If the limiting procedure is taken along a null direction instead, we get a notion of mass at null 
infinity known as the Bondi or Bondi-Sachs mass |59( 160) . This mass is time-dependent but 
never increasing, which means that it takes into account the energy that is carried away by 
gravitational waves (which is itself always positive). Furthermore, the positivity of the Bondi 
mass itself, and of the ADM mass, have been proved for physically reasonable space-times 
(PI and references therein). 


Fig. 4.1 illustrates the difference between the Bondi and ADM definitions. While the ADM 


mass is computed on the spatial (blue) slices, and thus measures the total (constant) energy 
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Figure 4.1: Schematic illustration of the different slices where the ADM and Bondi masses are 
computed: whereas the ADM mass is computed on a spatial slice (in blue), the Bondi 
mass is defined on a null slice (in red) and thus is sensitive to gravitational radiation (in 
yellow). From m- 


available in the space-time, the Bondi mass is computed on null slices (red cones) which previ¬ 
ously emitted radiation never reaches. Therefore, the Bondi energy is the energy remaining in 
the space-time at some retarded time, after the emission of gravitational radiation, and differs 
from the ADM mass precisely by the integral of such energy flux. All these three definitions 
of mass agree wherever their applicability regimes overlap (for instance, for a Schwarzschild 
black hole). 


Finally, another notion of energy-momentum exists in linearised theory. The Landau-Lifshitz 
pseudo-tensor m is essentially tI^) in Eq. ( |3.3.9[ ), i.e. it is an effective energy-momentum 
tensor generated by first-order metric perturbations. Despite not being unique nor gauge 
invariant, a well-defined gauge-invariant energy can be defined by a suitable integral at 
infinity. Its generalisation to higher dimensions [62] was actually the method employed in 
|54| for the first-order calculation. However, in what follows we shall focus exclusively on the 
Bondi formalism which is valid non-perturbatively. 


4.2 Bondi mass in D dimensions 

The Bondi formalism was only recently extended to higher dimensions [621 [S5] . The starting 
point is the choice of Bondi coordinates {T,f,x^), where f is a retarded time, f is a radial 
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(null) coordinate and are D — 2 angles, such that the metric can be written in Bondi form, 


I 


ds^ = grfdN + 2g^f.dTdr + 2g^^dx^df + gfjdx^dx 


The radial coordinate f is chosen to be an areal radius, i.e. such that 


(4.2.1) 


\9fj\ = "^^0-2 , 


(4.2.2) 


where ojd -2 is the volume element of the unit {D — 2)-sphere (not be be confused with the 
total volume il£)_ 2 ). In these coordinates, null infinity is located at f = oo. 

Asymptotic flatness is then defined by imposing an appropriate fall-off behaviour with out¬ 
going boundary conditions at null infinity, 

[fc-i-i] 


90 _ 
^2 


- ^ij + ^iJ + ^ 


f) 


ij 


k>0 


f.D/2+k-l 


(4.2.3) 


where ujfj is the metric on the unit (11 — 2)-sphere and k runs over all (non-negative) integers 
for D even and semi-integers for D odd. 

This, together with Einstein’s equations (in vacuum), is enough to fix the asymptotic be¬ 
haviour of all metric components. In particular, 


k<^ 


grf 


-1- E 


Ufc+i m{f,x^) 


-n r 2 


k=0 


‘-+k 


-b 


sD-3 


+ o(f-(^-i)) 


(4.2.4) 


The coefficients ak+i are unimportant but m{T, x^) defines the Bondi mass, 

D-2 I 


Mb{t) = 


m{f, p) dD,D -2 ■ 


16ttGi) Jbd-2 

Einstein’s equations guarantee its finiteness and yield an evolution equation, 

dM]^ 1 / 


(4.2.5) 


df 




(4.2.6) 


where the latin indices I are raised and lowered with gjj, was defined in Eq. (4.2.3) and 
the dot denotes a derivative with respect to retarded time f. 

In the remainder of this chapter we shall apply this mass-loss formula to our problem of shock 
wave collisions and derive an expression for the inelasticity, valid non-perturbatively. Before 


diving into the details, however, we will specialise Eq. (4.2.6) to axisymmetric space-times. 
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4.2.1 Bondi mass-loss formula in axisymmetric space-times 


In axisymmetric space-times one can choose one of the D — 2 angles to be the angle with the 
axis of symmetry, such that the remaining D — 3 angles lie on a transverse plane. So we split 


= ((9, 


and the metric components only depend on 9. Moreover, 


(4.2.7) 


f )j jdxMx-^ = d9‘^ + 1)^^ sin^ 9 dnjy_^ , 


(4.2.8) 


and the areal radius condition, Eq. (4.2.2), becomes 


+ ^00) (1 + ^,^0) 


D-3 


= 1 . 


(4.2.9) 


Thus can be eliminated in favour of In particular, asymptotically we have 


00- 


f) 


ee 


-(^-3)fl00. 


(4.2.10) 


Then the integration over the transverse angles (j)'^ in Eq. (4.2.6) can be made, and we get 
the following angular power flux 


dMs _ {D — 2){D — 3 )Dd-3 


fp 2 


1 2 


(4.2.11) 


dfdcos9 327rG£) r->-i-oo l 

where p = rsin0 is the radius on the transverse plane orthogonal to the axis of symmetry. 
We now have all we need to define the inelasticity of the collision in Bondi coordinates. 


4.2.2 Formula for the inelasticity of the collision 

The inelasticity of the collision, e, is the ratio between the energy radiated and the total 
initial energy. Since each of the colliding shock waves has energy E, and the radiated energy 
is the difference between the initial and final Bondi masses, 

Mb{t) 


e = 1 — lim 

f^oo zE 


(4.2.12) 


Recall that the energy E relates to the energy parameter k through 

SttGdE 


K = 


^D-3 


(4.2.13) 
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and in Section 3.3 we have chosen units in which k = 1. From Eq. (4.2.11) we then conclude 


that the inelasticity is given by the following time and angular integral 


dcosO 


€ = 


e(0) = 


dcos6 


dfs{f, oy 


(4.2.14) 


'-1 




where 


4 r^+oo 


D-4 • 

r/5 2 f) 


</></> 


(4.2.15) 


is essentially the generalisation of Bondi’s news function to higher dimensions. 

In the next section we shall study the relationship between de Bonder and Bondi coordinates, 


and rewrite Eqs. (4.2.14) and (4.2.15) in terms of quantities we know. We will see that Bondi 


coordinates actually asymptote to (a combination of) de Bonder coordinates, and that the 
news function is determined by the transverse metric perturbations hij only. 


4.3 Relationship between de Donder and Bondi coordinates 


In the de Bonder coordinates of Chapter the metric has the generic form 


ds‘^ = ds^i^^ + huudu^ + 2huidx^du + h^^dv^ + 2huvdudv + 2hmdx^dv + hijdx^dx^ . (4.3.1) 


First we change to Bondi-like coordinates (r, r, 0), where r = t — r is a retarded time, 
r = is the radius of a (B — l)-sphere and 9 is the angle with the z axis. The 


other angles (/>“ are ignorable so we omit them in this discussion. Eq. (4.3.1) now reads 


= dsyh^^+hrrdT"^+2hT-rdTdr+hrrdr‘^+2hT-QdTdd+hQQd9'^+2hr-edrd9+r‘^ sin^ 6hfp^dD,‘y_^ , 

(4.3.2) 
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where, from Eqs. (3.5.3), 


hrr = A + G + 2C , 

hrr = ^(1 — COS 6) + G{1 + cos 6) + 2C + {B + F) sin 6 , 

hrr = A{l-cos9f + G{l + cosef + {2C + H-{D-3)E)sEi^e + 

+2 (i?(l — cos6<) + F(1 + COS0)) sin0 , 

hre = r [(^ — G) sin 0 + (i? + F) cos 0] , (4.3.3) 

hee = [{A + G-2G)si-D^e + 2{B-F) sin 9 cos e + {H-{D-3) E)cos^ e] , 

hre = 'T [(^(1 ~ COS 9) — G{1 + cos 9)) sin 9 + (2(7 + H — {D — 3)E) sin 9 cos 9+ 

+ {B + F) cos 9 + {B — F) (sin^ 9 — cos^ 0)] , 

= E + H . 


Next we do the actual transformation to Bondi coordinates, 




(4.3.4) 


To make contact with Section 4.2, we denote the metric perturbations in these hatted 
coordinates by where the indices now run over all coordinates. The Bondi gauge is 
then defined by imposing = 1).^ = 0 and the areal radius condition, Eq. (4.2.2). 


Einding the vector that satisfies these three conditions is a demanding task. However, 
since we are only interested in the asymptotic behaviour at null infinity, we shall assume that 
decays sufficiently fast with some power of 1/f ~ 1/r, such that both coordinate systems 
match asymptoticall 3 Q Then the Bondi metric perturbations are given by 


and the solution (with = 0) is 

r = 


^+7(f,0), 


r = yy 


r = - 


h 


rr,e [Ke 7 ( f 0 ),^ . 

- / -^- + Pir,9) , 


2{D - 2) [ f2 


hee 




(4.3.5) 

(4.3.6) 

(4.3.7) 

(4.3.8) 


*We shall see later in Chapter|^that this is not true in D = 4. As discussed extensively by Payne 
contains Inf terms in D = 4. Here, however, we focus on D > 4, following |53| . 
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where f denotes primitivation with respect to f. 

7 (f, 0) and f3{f,8) are two arbitrary integration functions, constrained only by requiring the 


Bondi metric to be asymptotically flat, Eq. (4.2.3). In particular 

r 1 fh 


= fc** - ^ + 2(cot 9-d, 




(4.3.9) 


To have the required asymptotic decay, the contribution from /3 must vanish. Equating the 
differential operator acting on (3 to zero gives f3 = a{T) sin 0 . 


The same applies to the 7 contribution in Z) > 4. In iZ = 4, however, 7 remains arbitrary. 
This is the well known supertranslation freedom [421 , but we shall ignore it nonetheless 

(amounting to a choice of a particular supertranslation state). 


All that is left is to write Eq. (4.3.9) in terms of the scalar functions in Eqs. (4.3.3) (evaluated 


in Bondi coordiantes) and take a f derivative. The result, which depends on many of them. 


can be simplified by making use of the de Bonder gauge conditions, Eq. (3.3.8), which near 
null inhnity read 


, D-2 

lim f 2 

f^ + CXD 

. D-2 

lim f 2 

r—^-\-oo 


(1 — cos0)yl + sin 9B + ^^^{1 + cos 6)H 


4 

D-2 I 


(1 + COS 6)G + sinOF + ^^^{1 — cos 6)H 


lim f 2 

r^ + CXD 


(1 — cos 6)B + (1 + cos 9)F + sin 9 ( 2C — {D — 3)E — 


D-i 


H 


= 0 , 

= 0 , 

= 0 . 

(4.3.10) 


These, together with Eqs. (4.3.3), imply that 


. D-2 ■ . D-2 ■ 

lim f 2 Jirr = lim r 2 h^g = 0 

f^ + OO f^+OO 


lim f 2 


=2 hgg 


r^+oo 


= -{D-3){e + h] , 


which, inserted in Eq. (4.3.9), yields 


^^4. = E + H. 


(4.3.11) 

(4.3.12) 


(4.3.13) 


Therefore, the inelasticity is determined solely by the transverse metric perturbations. 

In > 4, we may drop the hats since (r, 9) —>■ (r, 9). However, in 0 = 4, the two coordinate 
systems do not match asymptotically. Nevertheless, we shall use the same expression, i.e. 
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work with de Bonder coordinates asymptotically. Later in Chapter we will see that this 
gives rise to logarithmically divergent terms in the second-order news function when r —)> oo. 
Then, following Payne [l5], we shall simply remove them. 


4.3.1 News function in de Donder coordinates 

In summary, if the transverse metric perturbation, in de Donder gauge, is decomposed into 
its trace and trace-free parts, 

hij = H5ij -\- ENij , (4.3.14) 


the news function defined in Eq. (4.2.15) reads 


'' 4 r—>.+00 


D-4 , 

rp 2 [E + H 


(4.3.15) 


where the dot denotes a derivative with respect to r and the relationship between (r, r, 6) 
and {u, V, p) is 


u = r-|-r(l — cos0), u = r -|- r(l -|-cos 0), p = rsm0. 


(4.3.16) 


Finally, the inelasticity, Eq. (4.2.14), becomes 


dcosO dcosO f , •, ^.n 

e= I -^e{e)= j dTS{T,ef. 


(4.3.17) 


'-1 


4.4 Integration domain 


The domain of integration in Eq. (4.3.17) is, in principle, r G ]—oo,-|-oo[ and 9 G [0,vr]. 
However, since no gravitational radiation is expected prior to the collision, one may ask if 


the news function in Eq. (4.3.15) is actually non-vanishing in that entire domain. This is 


particularly important if one wants to compute, and integrate, the news function numerically. 
Furthermore, since the transformation to Brinkmann coordinates adapted to the u-shock 
broke the Z 2 symmetry of the problem, one should check what effect, if any, does that have 
on the 9 domain. 
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4.4.1 Time domain 


Regarding the time domain, let us go back to Fig. |3.4[ The radiation signal at point V starts 
when it is hit by ray 1 and is expected to peak around the time it is hit by ray 2, which 
comes from the other side of the axis and has already crossed the caustic at p = 0. 


The trajectories of such rays are given by Eqs. (2.3.2) with u = 0 (and k = 1), 


V = <h(p) + u- 


p = ±p { 1 + u 


2p 


(4.4.1) 


where the + sign gives ray 1 and the — sign gives ray 2. 


Thus if the observation point V has coordinates (r, 9), the times ri and T 2 at which it is hit 
by those rays are the solutions of 


(2 4> -/>4>') =F 2pcot 5 e_ 

T-{p) = - - -Tcot-p + 0(4»), 

2 ± cot ^ 4>' 2 


r{p) = 


($'2 _Ap_ 2$$' 0 _ 


(4.4.2) 

(4.4.3) 


The asymptotic limit r —)• oo is achieved by p —)■ oo (we had already seen in Chapter that 
a far-away observer gets rays coming from the far-held region where the gravitational held is 
weak). Therefore it is clear that ri —)• —oo and T 2 —>• -|-oo for all D unless 6 = tt: an observer 
at the axis gets both rays at r = —oo va. D = 4 and at r = 0 in > 40 

The past light cone of an observer sitting on the axis at a hnite distance r is represented in 


Fig. 4.2, where it is manifest that the radiation signal will start as soon as it becomes tangent 


to the collision line v = <h(/o). Fig. 4.3 details these intersections for various points along the 
caustic. 


4.4.2 Angular domain 


With respect to the 0 domain, we note that the determinant of the transformation matrix 


given by Eqs. (2.3.2) is 


det 


d{u,v,p) 

d{u,v,p) 


l + ^ue{uW{p). 


(4.4.4) 


‘The discussion of this qualitative difference in |54| only applies to an observer exactly at the axis. In 
practice, we extract the radiation off the axis (even if only slightly), thus the integration limits to be considered 
are —oo to -foo, both in D = 4 and D > 4. 
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Figure 4.2: 3D space-time diagram (for the case D > 4) showing the caustic (in black) and three 
selected events, 1-3, on it. The past light cone of event 3 is drawn in red. Its intersection 
with the surface u = 0 (blue dashed parabola) is tangent to the collision surface u = 
0 A V = 4> (in green). The u-shock generators travelling along u that emerge from the 
intersection points will focus and converge at 3. The r axis is the world-line of an 
observer at hxed z, which will see no gravitational radiation before r = ri. The collision 
lines and the intersections with u = 0 of the past light cones of other points along the 

From \54^ . 


caustic are represented in the right panel of Fig. 


4.3 


Thus the transformation is well-defined for all h < 0, and for all h > 0 except at the axis 
p = 0. Actually, in the asymptotic limit f —>■ oo both coordinate systems coincide since 
p —>■ oo off the axis and hence —>■ 0. 

Furthermore, the radiation signal is expected to vanish for an observer at the axis since he 
sees no variation of the quadrupole moment. We shall see later that this is indeed true. 

Therefore we conclude that the domain of angular integration is 0 G ]0,'7r[. Given the Z 2 
symmetry, this should yield twice the integral over half the set (one of the hemispheres). 
Since our perturbative method consists of an expansion off 0 = vr, we anticipate that the 
limit 0 —)• 0 might be ill-defined at each order, in which case it would be appropriate to 
choose the —z hemisphere. 
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Figure 4.3: Collision surface, u = 0 A (green lines) and the intersection of past light cones 

(dashed blue lines) of various space-time points along the caustic with u = 0 for D = 4 
(left panel) and D = 5 (right panel). The past light cones are tangent to the collision 
surface. The latter case is qualitatively similar to any D > 5. The numbering 1-3 in 


the right panel corresponds to the space-time points with similar numbering in Fig. 4.2 
From 154^ . 









Chapter 5 


Angular dependence of the news 
function 


Symetrie. En ce qu’on voit d’une vue; 
fondee sur ce qu’il n’y a pas de raison de 
faire autrement. 

Blaise Pascal 
Pensees 


In the previous chapter we showed that the radiation signal seen by an infinitely far away 
observer, which when squared and integrated gives the energy radiated, is encoded in a sole 
function of retarded time r and angular coordinate 6. Furthermore, this so called news 
function depends exclusively on the transverse metric components hij, whose formal solution 
was given in integral form in Chapter 


In this chapter we shall simplify that solution and study the angular dependence of the 


news function. We start in Section 5.1 by showing that the metric possesses a conformal 
symmetry at each order in perturbation theory which, at null infinity, implies a factorisation 
of the angular dependence of the news function at each order. By actually computing that 
angular factor, we shall write the inelasticity’s angular distribution as a power series in 
sin^ 6. In particular, we show that a consistent truncation of this angular series at 0{n) 
requires knowledge of the metric perturbations up to 0(n + 1). This clarifies the meaning of 
perturbation theory in this problem: it allows for an angular expansion of the news function 
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off the collision axis. In Section 5.2 we simplify the integral solution and verify, explicitly, 
that it becomes effectively one-dimensional. 


5.1 A hidden symmetry 


D’Eath and Payne observed, in ZD = 4, that if one starts with a shock wave and performs a 
boost followed by an appropriate scaling of all coordinates, the end result is just an overall 
conformal factor |42) . This is easily understood in physical terms: since the only scale of the 
problem is the energy of the shock wave, the effect of the boost can be undone by rescaling 
the coordinates. In what follows we shall generalise this to ZD > 4. 


Recall the metric of one Aichelburg-Sexl shock wave in Rosen coordinates. 


ds^ = —dudv -|- ( 1 + 


ku6{u) 


dp‘^ + p‘^{l + 


Kue{u) 


2p 


4 >' 1 


(5.1.1) 


Now let L be the Lorentz transformation 


It (f> 


(n, V, X*) —> (e 


(5.1.2) 


and C the conformal scaling 


{u, V, X* 


C _^/3/- - 

—)• e V, X ). 


(5.1.3) 


Then, under the combined action of CL, 

(u,v,x ) —^ (e D-sPu,eo-3Pv,e (5.1.4) 

g^,(X) ^ eS^V(X'), (5.1.5) 

where X and X' denote the coordinates before and after the transformation respectively, i.e. 
X X'. The metric simply scales by an overall factor. 

Since this is a one parameter symmetry, one can find coordinates on (ZD — l)-dimensional 
sheets which are invariant under the transformation, and a normal coordinate parameterising 
inequivalent sheets. A suitable set of invariant coordinates on such sheets is {p, (?,</>“}, with 
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and (/)“ the angles on the transverse plane. The normal coordinate along the orbits of the 
symmetry is simply p which transforms as 


p = e 


D-3^ 


(5.1.7) 


Now consider the superposition of two shocks. In the boosted frame with Rosen coordinates, 
the action of CL is 

A; X) ^ X'). (5.1.8) 

Thus the perturbative expansion remains the same except for A —)• e“^^A. In the centre-of- 


mass frame with Brinkmann coordinates, Eq. (3.3.6) becomes 

r oo 

g^,{X) ^ g^,{X') = e-2^/'/ig(X') 

k=l 


(5.1.9) 


Quite remarkably, in the future of the collision the metric possesses a conformal symmetry 
at each order in perturbation theory. 

Note that the metric functions hlul in Eq. (5.1.9) are the same before and after the transfor¬ 
mation, the only change being the coordinates they are evaluated on. On the other hand, for 
a generic coordinate transformation, the metric transforms as a rank-2 tensor, 

dx^' dx^' 


gA^) = 


dx^ dx^ 




(5.1.10) 


Together with Eq. (5.1.9) this implies that 


/i(0(X') = , 


(5.1.11) 


where Nu and Aly are the number of u-indices and u-indices. Thus a given metric function 
evaluated on the transformed coordinates is the same function evaluated on the initial 
coordinates multiplied by an appropriate factor. 


In these coordinates, the invariants p, q read 


p = {v- ^{p))p 


D-A 


q = 


(5.1.12) 


while p transforms as 


V / — 7fQ/3 

p ^ p = e p. 


(5.1.13) 
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5.1.1 Reduction to two dimensions 


Since p is the only coordinate that transforms under the action of CL, Eq. (5.1.11) implies a 
separation of variables in the form 






(5.1.14) 


Since the angles (/>“ are ignorable, the problem becomes two dimensional at each order in 
perturbation theory. Let us see how. 


In Section 3.5 we defined scalar functions of {u, v, p) by factoring out the trivial dependence 

I - 1 (A:—1) 

of each metric component on Eqs. (3.5.3), and similarly for the sources T^u ■ Then 


Eq. (5.1.14) implies that, for each of those functions generically denoted by F{u,v,p) as 
before. 

For its respective source as defined in Eq. ( |3.5.5 ), since the d’Alembertian operator scales as 
p~‘^ under CL, 

s^’'\p,q) 


S^'^Hu,v,p) = 


p{D—3){2k-\-Nu—Nv)-\-‘2, 


(5.1.16) 


Thus our problem is effectively two-dimensional in (p, q) coordinates: at each order, q) 

is the solution of a differential equation (inherited from the wave equation) with a source 
s^^\p,q) and subject to initial conditions on = 0. This is an enormous computational ad¬ 
vantage. In Section [ b^ of Appendixj^we obtain the differential equation obeyed by f^^\p, q) 
and the reduced Green’s function G^{p,q;p',q'). However, there is another consequence of 
this symmetry which only becomes manifest when one goes to null inhnity. 


5.1.2 The CL symmetry at null infinity 

In terms of (r, r, 9) the new coordinates p and q read 

p = {t + r{l + cos9) — ^{rsm9)){rsm9)^~^ , (5.1.17) 

q = {t + r{l — cos6)){rsm9)~^^~^K (5.1.18) 

One can see they are not well-defined in the limit r —?• oo since 


p —)• oo 


(5.1.19) 
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However, it is reasonable to expect that the finite, non-trivial dependence of the news function 
in (r, 9) should be given in terms of combinations of p, q that remain finite and non-trivial at 
null infinity. An appropriate choice is 


p= < 


22^ + 21og9 , 

pq-1 

1 ) 
qD-3 


D = A 
D>A 


1 

q = qD-3 . 


(5.1.20) 


This coordinate transformation, {p,q) —)> {p,q), has a constant non-vanishing determinant, 

= D-3, (5.1. 

hence is well-defined everywhere. 

When r —)> oo, 


det 


d[p, q) 
d{p, q) 


p 2r(r,0) + C>(r-^) , (5.1.22) 

q —)■ (1 — cos6<)D^(sin0)“®^ X -|-O (r“^) , (5.1.23) 


where we have defined a new time coordinate 


r(T, 9 ) 


rx(l-cos0)-i + log(2^) , D = A 

1 D-4 

rx(l —COS0) 0-3 (sin 0) 0 - 3 ^ D>4 


(5.1.24) 


Recall that we had already factored out the p dependence in Eq. (5.1.15), leaving only a 
function of (p, q) to compute. But since r is the only surviving quantity in the asymptotic 
limit, we expect the news function to become effectively one-dimensional at null infinity: 
besides the trivial dependence on 9 coming from the known powers of r and p (remember 
that p = r sin0), it must be a function of r(r, 9) only. 


From Eqs. (4.3.15) and (5.1.15), 


d 


T('^)(r,0) = (sin0)'^-2^(^-3)x lim g) 

r—>-oo dr 1 


(5.1.25) 


where contains the relevant contribution from E^^'^ and as a function of {p,q). 

If the quantity inside brackets is to remain finite at null infinity, then necessarily 

+ , (5.1.26) 
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where ... denotes higher powers of q and a^^\p) some unknown function of p. 
Thus, after taking the limit r —?• oo, 

= -J— df ^ 

1 — cos V 1 — cos 9 dr dr 


(5.1.27) 


Using Eq. (5.1.24), and evaluating this expression at 0 = | noting that t (u f) = we 
conclude that 


9) = 


1 


1 — cos 9 
which proves our proposition. 


1 + cos 9 
1 — cos 9 


fc-i-U 


(r(T,0), 


TT 


(5.1.28) 


5.1.3 The meaning of perturbation theory 


The perturbative expansion of the metric, Eq. (3.3.6), implies an analogous series for the 


inelasticity’s angular distribution e{9) defined in Eq. (4.2.14), 

OO 

t(S) = ^ , 

where 


N=l 


N 

(0) = ^ / drT(^)(T,0)T(^+i-^)(r,0). 
k=i d 


(5.1.29) 


(5.1.30) 


Inserting the result of Eq. (5.1.28) and changing the integration variable from r to r, we 
conclude that 


e(^)(0) = 


1 


1 — cos ( 


3 /I , a\ ^-1 
1 + cos 9 


1 — cos 9 


:{N) A 


(5.1.31) 


Thus it suffices to compute the news function on the symmetry plane. 
The whole series reads 

3 OO / N-l 

1 + cos 9 


e(( 

N=1 ^ ^ N=1 ^ 

OO 

= y^^an{9)e 


— cos 9 


.(TV) 


TT 


n=0 


^(n+1) 

.2 


(5.1.32) 

(5.1.33) 


Observe that each individual 0^(0) is regular at 9 = tt, where we saw that perturbation 
theory should be valid, but does not obey an{9) = an{T^ — 9) since the transformation to 
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Brinkmann coordinates broke the Z 2 symmetry. However, e(0) should respect that symmetry. 
In particular, if it has a regular limit at the axis, it can be written as a power series in sin^ 9, 

00 


e{e) = . 


n=0 


Indeed, writing cos 9 = —y/l — sin^ 9, we see that near 9 = tt, 

an(^) ~ (sin0)^”. 


For example, the first three terms are given by 

« - (I) ■ 

ei = 


1 

J 2 1 

I 




€2 = 


128 




(5.1.34) 

(5.1.35) 

(5.1.36) 

(5.1.37) 

(5.1.38) 


Thus we conclude that a consistent truncation of the series in Eq. (5.1.34), i.e. the extraction 
of the coefficient e„, requires knowledge of and hence, from Eq. (5.1.30), of the 

metric up to 0{n + 1). This clarifies the meaning of perturbation theory in this problem: it 
allows for the extraction of successive coefficients e^, thus amounting to an angular expansion 
off the collision axis. 


It should be stressed that this result is a kinematical consequence of the CL symmetry only: 
the dynamical (wave) equations have not yet been used (except for the background solution 
of course). 


Einally, the inelasticity is then given by 


e 



dcos9 

2 


{sm9f^, 


E l ill 

( 2 n + l)!!^*" 

n=0 ^ ^ 

2 

60 + 361 + .... 


(5.1.39) 

(5.1.40) 

(5.1.41) 


In the next section we shall obtain simplified expressions for the integrals giving the asymp¬ 
totic metric functions that contribute to the news function. By working in Eourier space 
with respect to the retarded time r, we shall also confirm, explicitly, that the coordinate 


transformation of Eq. (5.1.24) factorises the angular dependence out of the integrals. 
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5.2 Asymptotic integral solution for the metric functions 


Back in Chapter]^ we wrote the formal solution for each metric function F{u, v, p) in integral 
form, Eq. ds:^ . Since we are ultimately interested in the news function S{t,9), we define 


the asymptotic waveform F{t,9), obtained from F{u,v,p), according to Eq. (4.3.15), 


F{t,9)= lim rp a —F{u,v,p). 
r—>-oo cIt 


(5.2.1) 


Actually, the Dirac delta in the Green’s function is better dealt with in Fourier space, so we 
define 

F{uj,9) = j dTF{T,9)e-^‘^U (5.2.2) 

This is equally useful to compute the inelasticity since, by the Parseval-Plancherel theorem, 


drF{r,9f = — du\F{u,9)\\ 


(5.2.3) 


The next steps are detailed in Appendix [C] due to their tedious and technical nature. 

In short, one needs to take the asymptotic limit and integrate in r to obtain the Fourier 
transform. Dropping primes on integration variables for ease of notation, we get 


F{u},9) = —i 2 / dpp 2 J d-4 (a;psin6*)5 ( uj 


+m'' 


1 + cos 9 1 — cos 9 




where Jiy{z) is the za-th order Bessel function of the first kind and 


;p , (5.2.4) 


S{x,y;p) = - du dve ™'^S{u,v,p). 


(5.2.5) 


Next, one finds it convenient to transform to a new frequency D, 


w —>• D = w 


— , ,£>-3 


sm \ ^ ^ 1 — cos 9 


together with a new waveform 


(5.2.6) 


9) = Vw'(D) F(a;(D), 9 ), 


(5.2.7) 


such that 


duj\F{uj,9)\^ = / dD\F{D, 


(5.2.8) 
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Then, using the CL symmetry, namely Eq. (5.1.16), one shows that, at each order k in 
perturbation theory. 




1 \ 2 /1 + COS 9 '' ^ ^ 




(5.2.9) 


1 — cos 9 J V1 “ CO® ^ / 
where the angular dependence is now completely factored out of the integral, and the (new) 
Fourier space waveform evaluated at 9 = which we shall abbreviate to 


.F(f2) = .F( 12,-1 , 


TT 


(5.2.10) 


is given by 

m) = - 


L»-3 


dRR^J^^^{2R)S{Q-\n-R). (5.2.11) 


This formulation in terms of the new frequency 12 will prove to be useful later on in Chap¬ 
ter for the evaluation of surface terms. However, for the time being, we can invert the 
transformation in Eq. (5.2.6) at 0 = ^, i.e. 


12^a; = 212B^, J^(12) ^ ^Q'(O) E (0(0)) 


(5.2.12) 


The net relationship between cu and cD is 

a) = w X (1 — cos0)o^(sin0)®^ , (5.2.13) 

which is equivalent, in real space, to a transformation of the time coordinate 

T —>■ f(T,9) = r X (1 — cos 9)~ (sin 9)~ . (5.2.14) 

Apart from a 0-dependent shift in Z2 = 4 (indeed an example of a supertranslation |53) i r is 


exactly the same as in (5.1.24). 


In the next chapter we shall specialise this asymptotic solution to the surface case and 
compute all the terms that contribute to the inelasticity and depend linearly on the initial 
data. This will allows us to extract the isotropic coefficient cq. 











Chapter 6 


Contribution from surface terms 


Shut up and calculate! 


David Mermin 


This famous quote is usually attributed to Richard Feynman, but is in fact due to David 
Mermin |64] . Indeed, a lot has been done and said in the previous chapters, but the physical 
quantities we wish to understand have not yet been computed. In this chapter we shall finally 
see some results by computing the surface terms, i.e. the linear contribution from the initial 
data to the news function and the inelasticity. This is the only term present at first order in 
perturbation theory, and is all we need to extract the isotropic coefficient cq. 


We begin in Section 6.1 with a review of the partial results of 


where the hrst-order 


waveforms for even D where computed numerically with a C++ code. These were later 
complemented by the results for odd D in [65], where an empirical fit formula was found 
for eo as a function of D, 



1 

D ■ 


( 6 . 0 . 1 ) 


In Section 6.2 we obtain all surface waveforms analytically and, as a corollary, prove that 


Eq. (6.0.1) is exact. Indeed, all contributions from surface terms to the inelasticity yield 
simple rational functions of D. However, for all but the isotropic term they are meaningless 
without the volume terms. Nevertheless, our results provide insight into an issue already 
encountered by D’Eath and Payne in D = 4 [421143j : starting at second order, both surface 


and volume waveforms have non-integrable tails at late times. In Section 6.3 we trace the 
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tE^^UKt _D = 4. e = 7r-0.001 



rpE^^UKr D = 6, 9 = tt - 0.001 




Figure 6.1: Waveforms in a rescaled r coordinate for even D: the signal starts at r = ri and peaks 


at T = T 2 , which are the retarded times corresponding to rays 1 and 2 of Fig. |3.4| From 

m- 


origin of these tails to the Green’s function and show that they are generic and may be present 
at all orders. D’Eath and Payne computed the news function at second order in 11 = 4 and 
concluded that the tails cancel when the surface and volume terms are summed, hence we 
expect a similar cancellation in H > 4. Again, many details are left to Appendix [E) 


6.1 Review of numerical results 


In 


a numerical method was set up, and coded in C++, to evaluate the surface integral, 


i.e. the first term in Eq. (3.4.3). At first order the trace of the transverse metric perturbation 


vanishes so there is only one scalar function, E, to compute. The asymptotic limit was not 
taken analytically. Instead, the waveform was computed with 9 = tt — 0.001 (it vanishes if 
0 = TT at finite r) for several (large) r, and the limit extracted numerically. 


The results are summarised in the hrst row of Eig. 6.1, which shows plots of the (rescaled) 


waveform for D = 4,6,8. As discussed in Section 4.4. H the signal always begins upon 
reception of ray 1 at ti and peaks at T 2 when the second optical ray arrives. The number of 
oscillations in the interval At = r2 — ti increases with D. 

At the time, technical difficulties prevented the numerical integration for odd D. These were 


later overcome and the picture completed in [65j . here shown in Eig. 6.2 


To get the inelasticity, another numerical integration was required. The results are sum- 
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Figure 6.2: Waveforms in a rescaled r coordinate for odd D: the signal starts at t = ri and peaks 


at r = r 2 , which are the retarded times corresponding to rays 1 and 2 of Fig. |3.4| From 

m- 


marised in Table 6.1 and Fig. 6.3 together with the apparent horizon bound, Eq. (1.3.1), 
for comparison. Interestingly, not only is cq < cah, but also both seem to converge when 
D —)• oo. Moreover, the numerical points fit nicely (within the numerical error < 0.1%) with 


a remarkably simple formula, Eq. (6.0.1). 


D 

4 

5 

6 

7 

8 

9 

10 

11 

CAH 

29.3 

33.5 

36.1 

37.9 

39.3 

40.4 

41.2 

41.9 

€0 

25.0 

30.0 

33.3 

35.7 

37.5 

38.9 

40.0 

40.9 


Table 6.1: Comparison between the apparent horizon bound, cah, and the first-order estimate, eg, 
obtained numerically. Note that eg < cah- 


We shall not give further details of the numerical methods and coding involved. A pedagogical 
presentation can be found in m- Instead, in the next section we will obtain the news function 
analytically and show that the fit formula for eg is exact. 
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AH bound &: First order estimate 



Figure 6.3: Plot of the apparent horizon bound, eAu (in blue), together with the first-order numerical 
estimate (red points). Observe the perfect fit with the formula for cq in Eq. 6.0. 1| (red 
line). From 


6.2 Analytical evaluation of surface terms 


The initial data on n = 0 were obtained in Chapter]^ In particular, from Eqs. (3.3.3) we see 
that they all have the generic form 


F(")(0,n,p) = /r(p)(n - ‘h(p))"0(n - $(/,)). 


( 6 . 2 . 1 ) 


However, the de Bonder gauge condition, Eq. (3.4.8) is yet to be imposed. In Appendix [P] 
we show that the initial data for the transverse metric components, i.e. the scalars E and 
H, are unaffected (actually, at first order, none of the components change), so we may move 
on. Details of the following calculation can be found in Appendix |E.1[ 


We start by feeding Eq. (6.2.1) to the simplified integral in Eourier space, Eq. (5.2.11), and 
compute 


Next, one finds it convenient to Eourier transform again, 


( 6 . 2 . 2 ) 


E{t) = — / dPLF{0)e 

I 


iflt 


(6.2.3) 
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Given the non-linear relationship between oj and n, this new t is not proportional to r (except 
in D = 4) but, for all practical purposes, gives an equivalent representation of the waveform 
(not to be confused with the Minkowski time). 


This was not necessary to compute the inelasticity. The integral could already be made in LI 


as in Eq. (5.2.8). The advantage is that J-{t) can be written in closed form, 

1 fc-i 




D-3 




1 d 
^'{R)dR 


<I>\R)-^R'^J^^j2R)fl,^\R)) , 

(6.2.4) 


where i? = <h ^(f). 

So, for each E^^\ one must compute the corresponding E via (6.2.4), which we denote 
by 

^ (6.2.5) 


Then, from Eqs. (4.3.15) and (5.1.30), the contribution to e(^^(|) {N = 1,2,3 for surface 
terms), is 


TTN {D-2){D-3) 


N 


E) A] - 


j dt [EE{k){t) + EH(k){t)\ [E^{N+i-k){t) + EH(.N+i-k){t)] , 

( 6 . 2 . 6 ) 


k=l 


where the bar denotes complex conjugation. 

Finally, we perform the time integration as an integral in the auxiliary variable R, 


j dtEjj(k) {t)E^(N+i-k){t) 



dR\^'{R)\EHik){<L{R))^E(^+^-^)iHR)), 


and similarly for the other contributions. 


(6.2.7) 


Thus we arrive at the striking conclusion that all surface integral contributions to the 
inelasticity are given by integrals over Bessel functions. 


Table 


6.2 


summarizes 


the contribution to e^^^(|) from each term in oc {E + H)^. Details 


of the relevant integrals can be found in Appendix E.2 


The first-order estimate for the inelasticity comes from N 
determines the isotropic term, 


1 which, from Eq. (5.1.36), 



1 1 

2~ D' 


( 6 . 2 . 8 ) 
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N 

Term 

contribution to e^^^(|) 

1 



2 

2-^ff(2) 

82 If d-3 

'^^\2 d) D-A 

3 

^ e;(2) 

2 ^” ^{2)Tfj(2) 

64(1 1) 

dJ {D+2){D+A) 

64(1 1)11-3 

l2 DJ D+2 

64(1 1) (^-3)" 

11^12 dJ {D-A){D-8) 


Table 6.2: Surface integral contributions to the news function series used to compute the inelasticity, 


Eq. (6.2.7). 


Thus we prove that the numerical fit was exact. Indeed, all the terms in Table 6.2 have been 
computed numerically with the same code used in [5l], and the results agree with a relative 
error of less that 10“^. 

For N = 2,3, observe that they are all rational functions of D, and all contain cq- Moreover, 
the minus sign at = 2 suggests that the series might be alternate, though this is of limited 
validity without the volume terms. The most striking observation, however, is that some 


terms have poles at Z) = 4 or D = 8. Indeed, in Appendix E.2 we show that some of the 
integrals do not converge for all D > 4. In the next section we will understand why. 


6.3 Late time tails 


Let us examine the late time behaviour of J-{t) in Eq. (6.2.4). The far future t —)> oo 
corresponds to Z? —)■ 0 since 


<1> ^(t) oc < 


6 2, D = 4 
, D > 4 


(6.3.1) 


For the surface terms in E and ZZ we get 


Eh( 2) (t) ~ R 


D-& 


(6.3.2) 


which means that IFu{ 2 ) grows exponentially with e* in ZZ = 4. Indeed, only for D > 8 does 
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it decay faster than which is the condition for integrability. This is in agreement with 
the hndings of D’Eath and Payne |43j in H = 4: both the surface and volume terms have 
non-integrable, exponentially growing tails at late times but their sum is well behaved and 
integrable. 

In Appendix |E.3| we trace the origin of these tails to the Green’s function and show that they 
are generic. In particular, we confirm that the volume terms of and have the same 
behaviour as their surface counterparts above. Therefore we expect a similar cancellation in 
D>4. 





Chapter 7 


The two-dimensional problem 


Odd as it may seem, I am my remembering 
self, and the experiencing self, who does my 
living, is like a stranger to me. 

Daniel Kahneman 
Thinking, Fast and Slow 


In the previous chapter, we managed to compute, analytically, all the surface terms at null 
infinity and their contribution to the inelasticity. This allowed us to extract the isotropic 
coefficient cq and thus obtain a first-order estimate for the inelasticity. At higher orders, we 
need to take into account the non-linear volume terms as well. These higher-dimensional 
integrals do not appear to be calculable analytically (actually, the same applies to the surface 
integrals if the asymptotic limit is not taken). Thus, numerical methods must be used. 


In this chapter we shall begin preparing the problem for numerical integration in a com¬ 
puter. This requires, amongst others, careful study of the integration region and choice of 
appropriate coordinates. 


We start in Section [7T| by using the CL symmetry to recast the general solution as a two- 
dimensional integral (one-dimensional for surface terms), which is an enormous computational 
advantage. In Section |7.2| we study the integration domain and identify the regions where 


singularities might be present. Then, in Section 7.3, we perform an hyperbolicity analysis and 


find characteristic coordinates which, after compactification, produce a finite two-dimensional 
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domain, presented in Carter-Penrose form in Section 7.4 Finally, in Section 7.5, we consider 
an observer at null infinity and further simplify the Green’s function. This is especially 
relevant for the second-order calculation of the next chapter. 


7.1 Reduction to two dimensions 


We start by recasting our dynamical problem in a two-dimensional form. From Eq. (3.5.6), 


we factor out the p dependence using Eqs. (5.1.15) and ( |5.1.16 ) (details in Appendix (B.3)) 
We get that the solution to each reduced metric function q) can be written as 


= J dq j dp'G^{p,q;p',q')s‘^^\p',q') 

where is the reduced Green’s function. 


(7.1.1) 


Giip,q-,p',q') = -J 


(7.1.2) 


and X* now reads 


l + y'^-{q-q'y^ ‘^){p-p'y - T(y)) 


X* - 

2y 


' 4>(y) , D = A 

IS 

III 

2 


Hy) . o>4 


(7.1.3) 

(7.1.4) 


Thus all the quantities that we might need, at any order in perturbation theory, can be 
computed numerically as a (at most) two-dimensional integral. Surface terms are a particular 


(and simpler) case. In Appendix B.3.1 we show that they are given by a one-dimensional 
integral with a structure similar to the Green’s function itself, 


where now 


xs = ®*(p' = 0, g = 0) = 


l + j/^-g(p-T(j/)) 

2y 


(7.1.6) 


The causal structure is encoded in the functions Im'^ixic) and in x* itself. In the next section 
we shall see how they determine the domain of integration. 
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7.2 Integration domain 

The integration domain is defined by two causality conditions. For an observation point 
V = {u,v,x^) in the future of the collision, the integration point V = must 

1. be in the future light cone of the collision; 

2. be in the past light cone of the observation point V. 

We shall study these two conditions separately. 


7.2.1 The future light cone of the collision 


Regarding the first one, recall from Chapter that, in Rosen coordinates, it is defined as 

h = 0Ai;>0 V h>0Ai; = 0. (7-2.1) 


In Brinkmann coordinates, from Eqs. (2.3.2), these become 
tt = 0 A u > ^{p) 




V 


u > 0 A V = ^{p) + 


(7.2.2) 


In the (p, q) plane, these two conditions define two important curves. The first is simply 
where the initial data have support, 

(7.2.3) 


p > 0 


A 


q = 0, 


while the second separates the flat region from the curved region, 

C“1 „ 


P(C) = ^(C) + 


D-3 


A g(C) = (C-i)C 


(7.2.4) 


where C £ [l)+c>o[ is a parameter. This is actually ray 1 of Chapter or the + sign in 

3y 

9(C) = (C + 1)C^-^ 


Eq. (4.4.1). The — sign, corresponding to ray 2, is similarly given by 

C + 1 


p{C) = 'I^(C) + 


D-3 


A 


(7.2.5) 


This curve is also important because the source is expected to be singular on it. Thus, extra 
care needs to be taken to integrate near this region. 


Einally, we conclude that the lower bounds for the integration variables in Eq. (7.1.1) are 

g'>0, p'>\ ^ . (7.2.6) 


—oo , 

D = A 

2 

D>A 

T)-4 ’ 
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7.2.2 The past light cone of the observation point 

For the second condition, we also begin the analysis in the more intuitive (u, v, p) coordinates. 


7.2.2.1 Three-dimensional case 


As mentioned in Chapter and explained in Appendix the original Green’s function 
G(ri, u, X*) has support on the light cone (y = 0) for even D but also inside it (y > 0) for 


odd D. From Eq. (3.5.8), these conditions are equivalent to — 1 < x* < 1 for even D and 


X* < 1 for odd D (indeed this is the domain where Im^ is non-vanishing). The region x* > 1 
is outside the light cone (y < 0) and thus causally disconnected from the event (u, x,x*). 

If we define curves (indeed parabolas) C±{p') by 

C±(p') = (p ± p')' -iu- u')iv - v') , (7.2.7) 

we can write the above conditions as 


C-{p') <044>x*-l<0, G+(p') > 0 X* + 1 > 0 , 


(7.2.8) 


In this way we can characterise the domain of the p' integration in Eq. (3.5.6). Both curves 
start at 

C±(0) = ~ {u — u'){y — v '), (7.2.9) 


which can be either positive or negative (or zero). Then C-\.{p') is monotonically increasing 
while C-{p') has a minimum at p' = p, 


C.{p) = -{u-u'){v-v')^-{Npf. 


(7.2.10) 


If this minimum is positive, the domain is empty. Thus {u — u'){v — v') >0 is an absolute 
condition for the event {u',v\x'^) to be causally connected with the event (u, x,x*). As this 
quantity increases, the domain is p' G [p — Ap, p-|-Ap]. Note that both ends correspond 
to X* = 1. The other condition, C+ > 0, begins to matter when {u — u'){v — v') = p^, 
i.e. when both curves start at the origin. After that, the domain is p' G [Ap —p, Ap + p], 
corresponding to x* = — 1 to x* = 1. 

Even when the only condition is C-{p') < 0, the location of G+(0) = 0 is important as it 
corresponds to a coordinate singularity at p' = 0: the angle (/>' between x* and x'* is not 
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defined if p' = 0, so neither is x* = cos cj)'. This is the first time that the p' integration starts 
from the origin, and it corresponds to the first ray that came from the other side of the axis. 
The Green’s function is expected to be singular there, and proper coordinates will need to 
be chosen to perform the integration near this region. 


7.2.2.2 Two-dimensional case 

In the two-dimensional reduced case, the picture is less intuitive but very analog. We define 
similar curves C±{y), 

C±{y) = (y ± - (y - qy^~‘^){{p - - p) , (7.2.11) 

such that, as before, 

C_(y) < 0 4^ X* - 1 < 0, C'+(y) > 0 44 X* + 1 > 0. (7.2.12) 


This time, both curves start at a non-negative value, (7^(0) > 0, and grow to infinity for 
large y. Moreover, C'+(y) > C-{y)\/y G Mg'. So, as before, the domain starts when C-{y) 
is tangent to the y axis. After that, the condition C-{y) < 0 gives a finite domain for the y 
integration. The other condition, C'+(y) > 0, comes into play also when C+{y) is tangent to 
the y axis. After that, the y domain is broken in two. 

The first condition defines the boundary of the light cone {p', q') of the event {p, q), whereas 
the second gives the location of the singularity of the reduced Green’s function G^{p, q‘,p', q'). 
They are the solutions to 

G±(y) = 0 A ^C±{y) = 0. (7.2.13) 

ay 

For each case there are two solutions, parameterised by y G M)]’ and labeled by n = ±1, 


where 


/ 

P 

/ 

q 


y^ ^ [p - T(y) - (1 ± y)Anip, y)] , 


1 


y 


D-2 


i±y ' 

An{p,q)_ 


A±{p, q) 


IT Vl + (2 + (D-4)p)(D-% 
{D - 2)q 


(7.2.14) 

(7.2.15) 


(7.2.16) 
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Figure 7.1: Characteristic curves in the {p, q) plane, showing the future light cone of the collision 
(in blue), the past/future light cones for two observers (in solid/clashed grey), and the 
singularities of the source (in red) and the Green’s function (in green) for observer 02- 
In = 4 (left diagram), the blue curve on the left goes to p = — oo, whereas in U > 4 
(right diagram) is starts at p = — 2/(ZI — 4). 


For the — solution, A± gives the two characteristics going through (p, q) which delimit the 
light cone. For the + solution (the curves where the Green’s function has a singularity), one 
can check that the one with A+ is inside the past light cone whereas the one with A_ is 
inside the future light cone of (p, q). 


In Fig. |7.1| we illustrate all these curves in the (p, q) plane for D = 4 and D = 5. The collision 
occurs at (p, q) = (0, 0) and the blue curves delimit its future light cone: the one to the right 
is where the initial data have support; the one to the left (ray 1) goes to — oo in D = 4 and to 


a constant in > 4, as in Eq. (7.2.6). Then we have two observers, Oi and O 2 , together with 
their past (solid grey) and future (dashed grey) light cones; their interiors are, respectively, 
to the left and to the right of those curves. Finally, the red curve corresponds to the second 
optical ray (ray 2), where we expect the sources to be singular, and the green curve is the 
location of the singularity of the Green’s function inside the past light cone of observer 02- 


It is possible to write these curves in a non-parametric form. For that, it is useful to define 
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yet a new set of coordinates (P, Q) through 




p , 

D = 4 

2 

”+0-4 ■ 

D>4 


Q = (D - 2)i2 + (D - 4)p)q . 


(7.2.17) 


Then we can solve for the parameter y, 

+ Q — n{D — 3) 


y = 


y/l T Q' i n{D — 3) 

Inserting back in the parametric equations, we find that they can be written as 


(7.2.18) 


C±n(P',Q') =C-n(P,Q), 


(7.2.19) 


where 

C±{P,Q) = < 


To recapitulate, the light cone (P', Q') of the event (P, Q) is defined by 

C±(P',Q')=C±(P,Q), (7.2.21) 

and the points {P', Q') where the Green’s function for the observation point (P, Q) is singular 
are given by 

C^{P\Q')=C^[P,Q). (7.2.22) 


P-1 


+ In 


\/l + Q i 1 


1 


P 


\/l + Q i 1 


2 J ^/TTQTI ’ 
^ ^/l + Q ±{D-3) 


D = 4 


D-2 


- (P - 4) 0-3 ^ D > 4 


(7.2.20) 


In the next section we will find characteristic coordinates and confirm that all these curves 
are indeed characteristics of the reduced differential operator. 


7.3 Characteristic coordinates 

For our initial problem, the coordinates (u, u, x*) were adapted to the initial data given on the 
characteristic surface u = 0. In the two-dimensional reduced version, it would be useful to 
find characteristic coordinates (^, rj) such that the principal part of the differential operator 
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contains only the mixed term d^d^i, and the initial data is on a characteristic line given by, 
say, constant 


In Appendix B.3 we obtain the differential operator acting on f^^\p,q), Eq. (B.3.2). The 


terms with highest derivatives come from 


- Idpd, + ((2 + {D- A)p)dp -{D- 2)qdgf + ... , (7.3.1) 

where we discard first derivatives. Transforming to {P,Q), 

dp^ dp + QP-^dQ, dg^ {D-2)iD-4)PdQ, (7.3.2) 


this becomes 


(D - AfP^dj, + 4Q{Q-{D- 2){D - 4)) 8^ - A{D - A)P{Q + D- 2)dpdQ + ... . (7.3.3) 


Now we make a coordinate transformation. 


(P,g)^(C(P,Q),7y(P,Q)) , 


(7.3.4) 


and the above operator takes the form 


fiZr^.Z^) 


dr] 8^ 8^ f dr] \ ^ ^ ^ 8^ ( 8^\^ 8^ 


8P 8P 8r]8^ drf' V^pJ 


where Zx = dqXj8pX and the characteristic polynomial is 

C{Z) = {D- 4)2p2 + 4g(Q -{D-2){D- 4))^^ - 4(P - 4)P(g + D-2)Z . 


(7.3.5) 


(7.3.6) 


The coordinates (.^,7?) are characteristics if C{Z^) = C{Zp) = 0, i.e. if they obey 


2g(g-(P-4)(P-2)) 8r] 
Q + {D — 2)(1 — -y/lT-j-^g) 8Q 
2g(g-(p-4)(p-2)) 8^ 
Q + {D — 2)(1 + y/TA-~Q) 8Q 


(P-4)P 

(P-4)P 


8r] 

~dP' 

8P ■ 


(7.3.7) 

(7.3.8) 


These clearly allow for a solution by separation of variables. Moreover, they reproduce the 
results of D’Eath and Payne |421135] in P = With a convenient choice of normalisation, 
a possible solution is 

e = C-(P,g), r/ = C+(P,g). (7.3.9) 

‘Note that (D — 4)P —>■ 2. 
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Thus we conclude that all the curves found in the previous section are indeed characteristics: 
for an event its light cone is delimited by 

^' = e, V' = V, (7.3.10) 


and the Green’s function is singular on 


e = ri, r?' = e, 


(7.3.11) 


respectively in the past and future of (^, r]). The integration constants were chosen such that 


r] = 0 goes through {p,q) = (0,0) (the left blue curve in Fig. 7.1) and ^ = 0 corresponds to 


the second optical ray (ray 2), previously obtained in parametric form in Eq. (7.2.5) (the red 


curve in Fig. 7.1). 


7.3.1 Compactified characteristic coordinates 

From the above results it is straightforward to conclude that the ranges of the characteristic 
coordinates in the future of the collision are 


^G]-oo,+oo[, r/G[0,+oo[, 


(7.3.12) 


with ^ < r] inside the light cone. This suggests introducing compactified coordinates, 


e = 




V^TTF’ x/TT^’ 

such that the integration domain becomes compact. 


(7.3.13) 


^£[-1,1], f)G[0,l]. 


(7.3.14) 


Then the volume integrals, Eq. (|7.1.1|), become 

ri 


[ di' r _ df,' 

J— 1 ^max{0,^'} 


d{p',q') 






where the Jacobian determinant is 


dip, q) 


dip, q) 


diP,Q) 


di^,v) 


8(^,1]) 


diP,Q) 


di^,p) 


diL fj) 




(D-3\ 

Q “2 L. 

P T^QT^ 

)-2 

)-3 

1 

1 



\D-A) 


y/l + Q 

(l-r)2) 

i (l-|2)l 


(7.3.15) 

(7.3.16) 

(7.3.17) 
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7.4 Conformal Carter-Penrose diagram 


With the compactified characteristic coordinates (^, 77 ) we can produce a conformal Carter- 
Penrose diagram, Fig. \7.2\ of the effective two-dimensional space-time, where the causal 
structure and the important curves we have been discussing become very clear. 


/+ 



Figure 7.2: Carter-Penrose diagram of the effective two-dimensional space-time. The curves and 
colours match those of Fig. 


7.1 


The collision occurs at = (—1,0) and its future light cone is shown in blue: ^ = —1 

contains the initial data and 77 = 0 is the first optical ray. The left boundary is the axis, 
p = 0, given here by ^ = r). Null infinity is at r) = 1 and future time-like infinity (/"’") 

is at = (1,1). 

The world line of a time-like observer is shown as a dotted black curve. When crossing the 
blue line it is hit by ray 1 and the radiation signal begins, peaking later when crossing the 
red line, ^ = 0, which corresponds to ray 2. 

The light cones at two instants, Oi and O 2 , are also shown, with solid and dashed lines for 






7.5. AN OBSERVER AT NULL INFINITY 


93 


the past and future branches, respectively. When these ‘reflect’ (effectively cross) at the axis 
the Green’s function becomes singular. The solid green line is the location of this singularity 
in the past light cone of O 2 . Note that this does not exist for Oi < 0) because ray 2 (the 
red line) is the hrst ray coming from the ‘other’ side of the axis. 

The radiation signal is extracted at , so the world line of a far-away observer is pushed to 
the right as f) —)• 1. That limit is the subject of the next section. 


7.5 An observer at null infinity 


Since we are ultimately interested in the extraction of radiation at null infinity, we shall see 
how that works in the two-dimensional picture. Indeed, the characteristic coordinates are 
well suited for such an asymptotic analysis. Taking the limit r —)• 00 with r and 9 fixed, we 
hnd that 




T(r,6») - 1 

1 


(£)_4)dV3 


D-I 


D = A 


t{t, 0) — 1 ] , T> > 4 


2{D-3) 1 , 

V ~ 0(r), 

{D - 9 


D-3 

(7.5.1) 

where r(T, 9) is the time coordinate in which the angular dependence of the news function 


factorises, Eq. (5.1.24). Thus we conclude that null infinity is when 77 — 7-00 for hxed and 
that the time integration at null infinity is equivalent to an integration in 


Regarding the Green’s function, observe that all the dependence on the observation point 
( 73 ,( 7 ) comes through x*. In {P,Q) coordinates, from Eq. (|7.1.3), we get, in = 4, 


X* = 


4-Q(P-P' + 21ny) , 4 + Q'(P - P'+ 2Iny)^ 

' y 


8y 


8 


(7.5.2) 


and, in P > 4, 


{D-2){D-A)-Q {D - 2){D - A) - Q\ 
2{D-2){D-4)y ^ 2{D-2){D-A) ^ 

, QP'/P . Q'P/P' L 

^2{D - 2){D - 4)7/^-3 2{P) - 2){D -if 


(7.5.3) 


In these coordinates it is manifest that x* is actually a three-dimensional quantity, depending 
only on Q, Q' and P — P' 01 P/P', respectively for D = A and D > A. 
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In the limit r —)• cx), recall that q ^ 0 and p —)• cx) with pg —)• 1. This implies 

P^oo, Q ^ {D-2){D-A). 


(7.5.4) 


Thus the term proportional to P diverges unless we scale the integration variable y. A natural 
choice is 

O 

D = A 


y = 


Q . 


(7.5.5) 


P' \ D-3 

y) y, d>4 


Then, asymptotically, 
1 

Xi, 

1 


^ 2 V <3 


A' 

2 


2 y 


D-3 


1 + 


Q'f{D-3) 2{D-3)y^-^' 


{D-2){D-4:) (d_4)A'b^ 


D = 4 


D>A 


(7.5.6) 


where 


A'D-3 = 


P'-P + 2 + - -21nQ 
O' 


P'\d^ ( 2{D-2)[D-3) 


P'- 2 ^ 


P' D-3 


D = A 


, 74>4 


(7.5.7) 


Q {D 2){D 4)y ("24 _ 4 )u-3^ 

Since x* —)■ a:*(A', Q'), if we factor out the trivial scaling of Eq. (|7.5.5), the Green’s function 


becomes effectively two-dimensional. 


GLiPN;p',Q') 



1 — (2 fc+A^ tj,—) 

Q'j 

) 

pi\ 


~P j 

) 


G^(A',Q'), 


D = 4 


(7.5.8) 


Gi(A',Q'), D>4 

Surface terms are a special case of the above, obtained by setting g' = p' = 0 in x*, but 


the same procedure follows. We scale y as in Eq. (7.5.5) but without Q' or P'. Thus, from 


Eq. (7.1.5) 




Ql-ik+N^-NpAO^^^^ 


D = A 


(7.5.9) 

In particular, this can be used to compute the asymptotic behaviour of the second-order 
source s^'^\p', q') when p' —)• oo. 
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7.5.1 The asymptotic metric functions 


From Eq. (7.5.8), since Q ^ r ^ in D = A and P 




D-3 


in D > 4, we conclude that the 


Green’s function decays at null infinity with a power of r (or, equivalently, of ry) equal to 

D-2 


- {D - 3)i2k + N^ - N,). 


Therefore, we define the asymptotic metric function as the finite limit 


(7.5.10) 


/(")(0 = lim /(")(f),0 X 

rj^oo 


Q 


l — [2k-\-Nu—Nv) 


D = A 


1 \ ^§E^-i2k+Nu-N.,) 


P 


, Zl>4 


1 ./max{0,f'} 


where 




d{p',q' 




X < 


^ \ 1 —(2/c+7Vii—A^ tj) 

Q' 


pi\%E^-0k+Nu-N,) 


D = 4 


(7.5.11) 


(7.5.12) 


, D>A 

From the point of view of numerical integration, this means that both the Green’s function and 
the source can be tabulated independently on a two-dimensional domain, prior to attempting 
the double integration. Indeed, there is a more natural choice of coordinates for the Green’s 
function. Let us define 


= C_(A',Q') = 
6rj = C+(A',Q') = 

together with their compactified versions, 

-5? 


<5^ = 


1 + {D - 
1 + {D - 4)^^ 

A- - 
07] — 


(7.5.13) 

(7.5.14) 

(7.5.15) 


’ -^1 + 6rf 

Observe the analogy with the definition of ry and ^ through C±{P'^ Q') (but note that, whereas 
P' > 0, A^ can be negative). For a given observation time the required region of (<5^, (Jry) 
is ^-dependent. What interests us, of course, is the full domain covered by all values of I,. In 
D = A this is simply 

(5^G[—1,0], (5ryG[—1,1], with 6 ^ < 57 ] . (7.5.16) 
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For D > A the picture is complicated by the fact that the denominator of Eqs. (7.5.14)-(7.5.13) 

1 

changes sign when ^ = —{D — A) u-s. One can show that it becomes 


(5^g[-1,0], 5?7G[^o, 1], with 5^ < 5fj, 


(7.5.17) 


plus an extra rectangle, 


5 ^ G [0,1] , (5i 7 G [—1, ^o] 


(7.5.18) 


In the next Chapter we will apply these results to the second-order calculation. In particular, 
we shall obtain, numerically, the second-order Green’s functions and sources for the metric 
functions and H^‘^\ 





Chapter 8 


The second-order calculation 


Nothing in the world can take the place of 
persistence. Talent will not; nothing is more 
common than unsuccessful men with talent. 
Genius will not; unrewarded genius is almost 
a proverb. Education will not; the world is 
full of educated derelicts. Persistence and 
determination alone are omnipotent. The 
slogan “press on” has solved and always will 
solve the problems of the human race. 

Calvin Coolidge 


In this chapter we conclude our journey with the first steps towards the computation of the 
second-order news function, and its contribution to the inelasticity’s angular distribu¬ 
tion, e( 2 )(|). 


In terms of the asymptotic functions defined in Eq. (|7.5.11|), the news function reads 


{D-2){D-3) d 


4 

so our task is to compute the second-order metric functions and 


( 8 . 0 . 1 ) 


This would allows us to finally extract the second coefficient in the inelasticity’s angular 
series, 62 , and thus obtain an improved estimate for the inelasticity of the collision. At the 
time of writing, we were not able to finish this demanding task. The results here presented 
are only partial but, nevertheless, constitute an important step towards the final goal. 
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The numerical implementation of the strategy described in this thesis was done in a C++ code 
by Marco Sampaio, building on the earlier code used in [5ll|65] for the hrst order calculation. 
Therefore, we will not attempt to give a thorough description of the coding involvedj^ 

We begin in Section |8.1| by computing and tabulating the reduced Green’s functions in their 
full domains. Then, in Section [8.2[ we do the same for the sources. Finally, in Section 
we comment on the outstanding two-dimensional integration, which we hope to perform in 
the near future. 


8.3 


8.1 The second-order Green’s functions 


We begin by tabulating the two reduced Green’s functions (m = 0 and m = 2) 


( 8 . 1 . 1 ) 


with X* as in Eq. (7.5.6) and A' as in Eq. (7.5.7). We shall do this in the compact {6^,6r]) 
coordinates. 


Before doing so, though, we must ensure that they are finite within the whole domain, 
otherwise numerical errors will lead to loss of precision when interpolating the functions near 
the singular points. In Appendix we show that the Green’s function remains finite as 
r/' —)• oo, but on the axis rj' = it vanishes for D = A and diverges when D > A (with a 
power law). Indeed, the strongest power occurs at the intersection of the axis with the light 
cone, In the new coordinates this corresponds to 5r] = 5^ and hence, to make the 

Green’s function regular, we multiply it by an appropriate power of ((5?) — <5^), which can be 


read from Eq. (E.2.23). 


Regarding the green line (ry' = ^), which now occurs at 5^ = 0, we have found empirically 
that the singularity there is at most logarithmic^ To regularise, we choose to multiply by 6fj 
(indeed this makes the function go to zero at the green line). 

‘The interested reader may find a pedagogical discussion on the evaluation of one-dimensional integrals 
in m- Note that these make up not only the surface terms (and hence the second-order source) but also the 
second-order Green’s function for the volume terms. 

Mhis could probably be proved analytically in a manner similar to the asymptotic expansions of 
Appendix [F| but we did not find it necessary to do so. 
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Finally, to compensate the power of Q' that is now part of (Eq. ( |7.5.12 )) in Z? = 4, and 
which would make the Green’s function grow close to the initial data ((i.^ —>■ —1), we multiply 
by (1 + (51)1 

So, putting everything together, we are going to tabulate numerically the functions 

Sfj) = Gl{A', QO X 5fi {6fj - (1 + 5^^ , (8.1.2) 


where 


7i = < 


7 

~2 

3- 


2{D - 3) 


, D = A 
, D>4 


12= { 


D = A 
D>4 


(8.1.3) 


In Figs. 8.1 and 8.2 we plot the numerical results for Qg and G 2 for D = 4,6 and D = 5,7 


respectively. In the horizontal axis we have 5fj and in the vertical axis 5^. The domain is as 
explained before at the end of Section |7.5.1[ As usual, the green line depicts the location of 
the singularity, which has however been regularised. 

In these plots we have also suppressed the late-time growth identified in Chapter In = 4 
we chose a combination which goes as close to {5^, 6fj) = (—1, —1), 


exp 


1 


- 5 ? 


+ 


Sfj 


(8.1.4) 


2 y 1-^,5^ 1 + Sfj^ 

Strictly speaking, this growth is only expected in Qq, but we have used it for G 2 as well to 
smooth out numerical errors very close to (-1,-1). 

For 4 < H < 8 this growth manifests itself close to Sr] = ^g. Hence we multiply Gg (only) by 

, 8-D 


((5f) - 4 ) 


(8.1.5) 


We have not included these factors in the definition of Gm above because later when calling 
the function for the two-dimensional integration we shall simply divide the tabulated function 


by them, thus restoring the definition of Eq. (8.1.2) 


8.2 The second-order sources 

We have not yet written expressions for the second-order sources. Recall from Chapter 
that the source for the wave equation obeyed by metric perturbations is an effective energy- 
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Go {D = 6) 


G2 {D = 6 ) 



Go {D = 8) G 2 {D = 8) 



Figure 8.1: Second-order Green’s functions Go (left column) and G 2 (right column), for D = 4,6,8, 


including the extra factors of Eqs. (8.1.4) and (8.1.5). The domain is explained in 


Section 7.5.1 The green line shows the location of the singularity (which has been 


regularised). 
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Go {D = 5) G2 {D = 5) 



Go {D = 7 ) G2 {D = 7) 



Go {D = 9) G2 {D = 9) 



Figure 8.2: Second-order Green’s functions Go (left column) and G 2 (right column), for D = 5,7,9, 


including the extra factors of Eqs. (8.1.4) and (8.1.5). The domain is explained in 


Section 7.5.1 The green line shows the location of the singularity (which has been 


regularised). 
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momentum tensor built out of lower-order perturbations, Eq. (3.3.9). At second order, we 
have 




^ajS^iTU 


,(i) 


-h 


(l)a 

u,l3 


^ At,a ,At^a/3,it ' 

( 8 . 2 . 1 ) 


Using Eqs. (3.5.3), we can then read the scalar sources for the transverse metric functions 


and H^‘^\ as defined in Eq. (3.5.4), 


tI,^\u,v,p) = + 2{D - 3)E^E, - 2EEp, 


+ 


I 4(D-3) p2 2{D-3)^ pp, 2{D-3) pp -l-AF -\- 
^ ^ (D-2)p^^^P D-2 ^-^,pp + ^^.vv + 


+{D - 4) [2E„E, - 

rj;^(n,n,p) = [^EE, + E,pE, + EE 


( 8 . 2 . 2 ) 


+ 


+(E-3) 


p - ,u i I - ,pv 

4(JI-3) c,2 I D /TP \2 D^-QD-He 


- 


zr TP D —4 

(D-2)p D-2^^^PP 


EE, 


- 2{D - 3)E„E, 


(8.2.3) 


Finally, their two-dimensional form in {p,q) coordinates, t^^\p,q) and P^{p,q), is obtained 


by applying the definitions in Eqs. (5.1.15) and (5.1.16), and noting that the derivatives 


transform as in Eq. (B.3.1). 


The computation of these sources starts with the reduced functions a(^)(p, q), b^^\p,q), 
e^^\p,q) and their derivatives, all of which are given by one-dimensional integrals as in 
Eq. ( 7.1.5| ). The initial data, i.e. the numbers ag^\ b'^\ e^\ are obtained from Eqs. (3.3.3) 
and from the transformation to de Bonder gauge in Appendix (which only affects ). 
Then, we build t^^\p,q) and t^^\p,q). 

As for the Green’s functions, we need to ensure that the tabulated functions are regular and 
smooth. We choose to include also the factor so we must consider its behaviour as 

well. 

Regarding the sources themselves, we found, in accordance with the asymptotical analysis of 


Section 7.5 (in particular Eq. (7.5.9)) that when r/ —)• oo both grow with 


(1 -^) 




(8.2.4) 


At the axis, fj = ^, we found numerically that both sources go to zero. 
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At the red line (^ = 0), we know from previous experience that the first order metric functions 
are square integrable, so the singularity (which is only present in even D) is regularised if we 
multiply by 

Regarding from Eq. 

—>■ ±1. In the first case, we have 

~ r/'~ (1, Q'-{D-2){D-A)^r]'~^ , (8.2.5) 

so we see that, at leading order, goes to zero as 

(1 _y^')i+i(^-4) , (8.2.6) 


(7.3.17) we expect some singular behaviour when r)' —?■ 1 or 


Combining with the divergent behaviour of the source, Eq. (8.2.4) above, and given that the 
Green’s function is finite, we conclude that the full integrand goes with 




(8.2.7) 


In H = 4 this gives rise to the logarithmic terms discussed by Payne [45], which we shall 
subtract and discuss in the next section. For the moment, we just multiply by (1 — fj') to 
make the tabulated function regular for all D. 

( ^ 2\ i 

In the second case, we found that multiplying by I 1 — I makes the result regular. 

So, putting everything together, we shall tabulate numerically the functions 




E,H\ 


{p\ g') X e' 1 - (1 - v'){v' - (1 + ■ (8.2.8) 


In Figs. 8.3 and 8.4 we plot 5o and S 2 , for H = 4, 6, 8 and D = 5,7,9 respectively, on the 


conformal diagram of Fig. 7.2, Note that the original singularity of the source at = 0 has 
been regularised, so the red line is shown for reference only. 


The apparent growth near the intersection of the axis with the red line is due to the extra 
factor that came from the Green’s function. However, in the next section we will see that 
their product is smooth and finite. 
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So {D = 4 ) 



So {D = 6 ) 


52 {D = 4 ) 


52 {D = 6 ) 





5o {D = 8) 


52 {D = 8 ) 




Figure 8.3: Second-order sources 5o (left column) and S 2 (right column), for D = 4,6,8. The 
axes match the conformal diagram of Fig. |7.2| The red line shows the location of the 
singularity (which has been regularised). 
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So {D = 5 ) 


52 [D = 5 ) 




So {D = 7 ) 


52 {D = 7) 




So {D = 9 ) 


52 {D = 9 ) 




Figure 8.4: Second-order sources 5o (left column) and S 2 (right column), for D = 5,7,9. The 
axes match the conformal diagram of Fig. |7.2| The red line shows the location of the 
singularity (which has been regularised). 
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8.3 The two-dimensional integration 


Comparing Qrn x Sm with the original integral, we see that now we have 
/( 2 )(^^) = di' C _ dfl'jgm{6C,6v)Sm{i',fl'), 

J— 1 ^max{0,^^} 

where 




1 


1 1 1 f T]' -g 


o' C' \ T'l 


i + i' 


72 


(8.3.1) 


(8.3.2) 


contains all the singularities of the original integrand: either they are integrable, or the 
principal value must be taken in accordance with the finite part prescription. The only 
major issue left is the logarithmic term in H = 4. As seen in the previous section, the limit 


hm (1 - f,')JG^{6i,6r,)Sm{i\f,') = C{i') 


(8.3.3) 


is zero in H > 4 but finite in H = 4. This gives rise to logarithmically divergent terms when 
r/' —7- 1. Payne [U] showed that these can be removed by an explicit transformation to Bondi 
gauge, so we can simply subtract them, i.e. in 11 = 4 we write 

dfj' 


I -/ — I ^max{0,^^} ^ ^ 




+ di' c{i') iog(i - i') - iog(i - r)) di' c{i') I, 


(8.3.4) 


and we discard the last term. 


In Figs. 8.5 and 8.6 we plot the product Gm x Sm, for H = 4, 5, 6, on the conformal diagram 


of Fig. 7.2, In Fig. 8.5 the observation time r is between ri and r 2 , i.e. before the second 


optical ray, whereas in Fig. 8.6 we have t > T 2 - 


These functions are smooth in the whole domain, so we expect this two-dimensional integra¬ 
tion to be performed in the near future. The major outstanding challenge is the achievement 
of the required numerical precision. 
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Figure 8.5: Product of the second-order Green’s functions with the second-order sources, Qo x So 
(left column) and Q 2 x ^2 (right column), for D = 4,5,6, for an observation time r 
between ti and r 2 , i.e. before the second optical ray. The axes match the conformal 


diagram of Fig. 7.2 The green and red lines show the location of the singularities (which 


have been regularised). 
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GoSo (D = 4,i = 0.5) 



GoSo (D = 5,i = 0.5) 



G 2 S 2 {D = 4, C = 0.5) 



G 2 S 2 {D = 5,i = 0.5) 



GqSo {D = 6,i = 0.5) 


G 2 S 2 {D = 6, = 0.5) 




Figure 8.6: Product of the second-order Green’s functions with the second-order sources, Qo x So 
(left column) and Q 2 x 52 (right column), for D = 4,5,6, for an observation time r 
greater than T 2 , i.e. after the second optical ray. The axes match the conformal diagram 


of Fig. 7.2 The green and red lines show the location of the singularities (which have 


been regularised). 









Chapter 9 


Conclusion and 


outlook 


Perfection is achieved, not when there is 
nothing more to add, but when there is 
nothing left to take away. 

Antoine de Saint-Exupery 

In this final chapter we conclude with some remarks and an outlook on open questions. 

We believe this work constitutes a significant contribution to the understanding of gravita¬ 
tional shock wave collisions, which have been studied by many in the past few decades, with 
varied motivations. 

It is quite remarkable that quantitative information from a non-perturbative process like 
black hole formation can be extracted using perturbation theory. The CL symmetry, though 
not an exact symmetry of the space-time, seems to play a major role to this end. The 
understanding of how it justifies perturbation theory, i.e. the correspondence between the 
perturbative series of the metric and an angular expansion of the news function is one of the 
major results of this work. 

Extra spatial dimensions may not exist, but the simple formula we obtained for the isotropic 
contribution, cq, illustrates the advantage of considering the number of space-time dimensions, 
D, as a parameter in general relativity. The known four-dimensional result, cq = 0.25, would 
never look so understandable and appealing if the D-dimensional computation had not been 
performed. 
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The results of Table |6.2| for the second-order surface terms suggest that higher-order terms 
may exhibit similarly simple patterns, with increasing powers of 1/D. If this is the case, it 
would be interesting to relate our angular series with the large D limit of general relativity 


These remarkable analytical results were made possible by formulating the problem in Fourier 
space. It is tempting to think that an analytical solution might exist for the volume terms 
as well, so this is certainly an avenue worth pursuing. Indeed, given all the symmetries and 
simplifications that our work has revealed, perhaps there exists a more sensible approach to 
the problem which allows for an analytical solution. 

The second-order calculation has proved to be very technical and challenging, both analyti¬ 
cally and numerically. We expect the first angular correction, ei, to be obtained very soon, 
shedding some light on the issues above. 

Moreover, we hope our work may as well be the foundation stone for further studies. The most 
obvious candidate is the third (and higher) order contribution. From the numerical point of 
view, there is an increase in computer effort since three-dimensional Green’s functions would 
have to be mapped, but the integrals from which they arise are essentially the same. 

As argued in Chapter the fact that the initial data are exact at second-order, and the good 
agreement of the four-dimensional result with numerical relativity, suggest that the second- 
order estimate for the inelasticity might be accurate enough. The third-order calculation 
could not only validate (or not) this statement, but also provide a better understanding of 
the angular distribution of the radiation. 

Perhaps more interesting, though, would be to apply this method to different problems. 
Charged solutions have been considered in [53] , but the results seem to indicate a breakdown 
of perturbation theory. The introduction of an impact parameter, however, is probably 
treatable. Scattering with large impact parameter has been recently addressed in m, but 
otherwise the best estimates come from apparent horizon bounds |39j . It would be particularly 
interesting to see the threshold of black hole formation as the critical impact parameter is 
exceeded. 


Appendix A 


The infinite boost 


This appendix contains details on the inhnite boost of the Schwarzschild solution that yield 
the Aichelburg-Sexl shock wave. It supplements Chapter 


When taking the limit u —)> 1 in Eq. (2.1.6), we expect a distribution similar to the delta 


function. Therefore, it is natural to take that limit on an integral of the distribution. Let us 
dehne 


Fy{t,z,p) = dzj {z'- vtY + p- 

J —oo *- 

rZy{z,t) 


^ r dz' (z«+1)- 

H J —oo 


D-3 

2 


with Zy(z^t) = — vt)p To take the limit, we first note that 


(A.O.l) 


lim Z.u{z, t) = sign{z — t) x oo . 

v—>-l 


(A.0.2) 


1 


so 

lirn Fy(t, z, p) D—i / d-s D—4 ” -n /" D—3\ ’ 

^ J-oo (z'2 + 1)— ^ p r(^) 

where T is the gamma function. 


«(2-*) =r(^) 

-f;— i—V TT-77=;—TV 


(A.0.3) 


Finally, in Eq. (2.1.6), 


r .^-2 

lim 4—- 

v^i D — 3 


Yz = 


WttGdE 


lim 


dFy SttGdE 


{D — 3)^}d—2 r"—>-1 dz fl £)_3 


^{p)6{t- z) 


(A.0.4) 
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Appendix B 


Green’s functions 


This appendix contains details on the Green’s function for the d’Alembertian operator and 
on the general integral solution to the wave equation with a source with characteristic initial 
data. It mainly supplements Chapterand ChapterIn Section B.I we present the Green’s 
function for the D-dimensional wave equation and show how to use it to solve a problem 


with characteristic initial data. In Section B.2 we perform the dimensional reduction to 


three dimensions and, in Section B.3 to two dimensions. Finally, in Section B.4 we present 
expressions for the functions Im’"' and their derivatives. 


B.I Green’s function for the d’Alembertian operator 


The d’Alembertian operator □ is defined, in cartesian coordinates (t,z,x*), as 

□ = = -df + d‘^ + 0*5*, (B.1.1) 

where r]^'^ is the inverse Minkowski metric. 

Its Green’s function G{t,z,x^) is dehned as the solution of 

□G(t, z, X*) = 6{t)6{z)6{x^). (B.1.2) 


Observe that this is not a coordinate-invariant definition since the delta function 5{x^) 
transforms with the Jacobian determinant of a coordinate transformation, 




dx‘ 


ja 


dx^ 


5{x 




(B.1.3) 
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In particular, for (t, z) —)■ {u, v) = {t — z,t + z), 

6{t)6{z) = 26{u)6{v). 


(B.1.4) 


The Green’s function in cartesian coordinates, i.e. the solution to Eq. (B.1.1), is 


G{t,z,x^) = - ^ 

27r 2 


2 ) {t^ — — x''Xi). 


(B.1.5) 


Therefore, the solution to = (5(x^) in our coordinate system, x^ = {n, u,x*}. 


G{u,v,x^) = - X = - 

47r 2 


rj^yx^x^' . 


(B.1.6) 


For odd -D, we have a fractional derivative of the delta function. This is to be understood in 
Fourier space, or formally as follows. For negative orders (n > 0), 

5^-^\x) = ^^x^-^e{x), (B.1.7) 

where T is the gamma function. Then, acting with derivatives, we can extend this recursively 
to positive fractional orders. For instance, 

5i^{x) = ^5i-^){x) = -^^ (^-^-5{x)\ . (B.1.8) 

ax \/AA V J 

This also shows that, while for even D the Green’s function has support on the light cone 
(x = 0)5 for odd D it has support also inside it (y > 0). An odd number of space-time 
dimensions behaves as a dispersive medium. 


B.1.1 Solution with characteristic initial data 

Given an equation □F(u, u,x*) = S{u,v,x’') it is straightforward to show that a solution is 

just by using the definition of G{x^). But what is the solution to □T(u, u,x*) = 0 given 
T(0, u,x*), i.e. with characteristic initial data? Theorem 6.3.1 of |55j says that this is 
equivalent to a source 

S{u, V, X*) = 4(5(u)d^T(0, v, x*). (B.1.10) 

For a motivation of this result in Fourier space, see the appendix of |54j . 
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B.2 Reduction to three dimensions 


In this section we specialise to the case of spherical symmetry in the {D — 2)-dimensional 
transverse space spanned by the coordinates x*. We choose spherical coordinates with radius 
p = -s/x^Xi and D — 3 angles (/>“, and define a basis of tensors in that space, 

dij, = Aij = 6ij-{D-2)rirj, (B.2.1) 

with ranks m = {0,1,2} respectively. Let XmiXi) be one of such tensors. Our problem 
becomes to solve the equation 

□ (F(ri, n, p)Xm{Ti)) = S{u, v, p)XmiTi). (B.2.2) 


First, observe that we can write 

A = d% = p-^^-^^dpip^-^dp) + p-^Asn-3 , 


(B.2.3) 


where Agos is the Laplace-Beltrami operator (or the spherical Laplacian) on the {D — 3)- 
sphere. Indeed, our tensors Xm are eigenfunctions of this operator. 


AgD-sXm = -m{m + D - i)Xm , 


(B.2.4) 


thus the equation for F{u,v,p) implied by Eq. (B.2.2) is 


[-Idudv + dp + {D - 3)p ^dp- m{m + D - 4:)p F{u,v, p) = S{u,v, p). (B.2.5) 


Instead of trying to find a Green’s function for this operator directly, we can use the integral 


solution of Eq. (B.1.9). Contraction with the appropriate tensor on both sides yields 


F{u, V,p)= I d^x' G{x - X%m ( ) S{u\ u', p') 


n / 

X x, 


(B.2.6) 


where Am{z) = {l, z,{D — 3) ^ [{D — 2)z‘^ “ l)}) respectively for m = {0,1, 2}. 

Then again we choose spherical coordinates inside the integral, with the angles aligned 
such that one of them is the angle between x* and x'*. Then we integrate the remaining D —4 
to get 

^ \ [ du'dv'dp' p'^~^S{u',v', p') 


F{u,v,p) = - 


D-2 

271 2 


du'dv'dp' p'^ ^S{u',v',p') I dx{l — x‘^) 2 ®+”*Am(x)(5( 2 )(x). 

(B.2.7) 
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Using the scaling properties of the delta function, and defining x* through 


X = {u — u'){v — v') — — 2pp'x) 

= 2pp{x - X*), 


(B.2.8) 

(B.2.9) 


this becomes 


F{u,v,p)= / du / dv' / dp'Gm{u — u ,v — v'] p, p)S{u ,v', p'). 


(B.2.10) 


where 


Gm{u-u',V -v'-,p,p') = 


(B.2.11) 


is the new Green’s function and we have defined one-dimensional functions Im"' through 
/^’'"(x*) = ^^0^2 [ dxAmix){l-x^)^~U(^~''){x-x^), 

{27r)— J-i 


Od- 


, N 0-2 

( 27 r) 2 


^ A-W dxa™(l-x2)V+™5(V-«)(x-x*), (B.2.12) 


1-1 


where = {1 ,-{D - 3), {D - 1){D - 3)}. 


B.3 Reduction to two dimensions 


In Chapter we show that the problem becomes effectively two-dimensional order by order. 
In this section we obtain the differential operator acting on f^^\p,q) and its associated 
Green’s function G^{p,q;p'q,q'). 


From the definition of p, q in Eq. (5.1.12), we see that 


du-^ P ‘^^dq , dy ^ p^ '^dp , dp^ dp - —— -qdg + — —dp . (B.3.1) 

r y y p -i pi' 


Inserting this in Eq. (^B.2.5), together with Eqs. (5.1.15) and (5.1.16), we conclude that the 


operator acting on f^^\p,q) is 


-4dpdq + ((2 + {D- I)p)dp -{D- 2)qdq - {D - 3){2k + Ny + Ny) + D - 4) x 
X ((2 + iD- 4)p)dp -{D- 2)qdq - {D - 3){2k + + Ny)) - m{m + D - 4). 


(B.3.2) 
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As in the previous section, its Green’s function can be found from the integral solution for 


F{u,v,p), Eq. (B.2.10). We factor out the p dependence using Eqs. (5.1.15) and (5.1.16), 


and change the integration variable p' —)■ py. We get that the solution to each reduced metric 
function {p, q) can be written as 


f^^\p,q) = j dq' j dp'Gt{p,q;p',q')s'^^\p',q), 
where is the reduced Green’s function, 

Gi{p,q-,p',q') = r 


(B.3.3) 


(B.3.4) 


and now x* reads 


i + y^-(g-g'y^ ^)(p-p'y ^^-^(y)) 


Xi, - 

2y 


' ‘h(y) , D = 4 

III 

2 


m o>4 


(B.3.5) 

(B.3.6) 


B.3.1 Surface terms 

Surface terms are a particular case of this general solution. Eor those, the source is as in 
Eq. (E.1.1), where fQ^\p) oc p-^(-D-2)-(Af„-Af„)(D-3)^ 


s^^Hp,g) = 2k\f^^\l)6{q)6‘^ ^\p). 


(B.3.7) 


Inserting this in Eq. (B.3.3), the q integration is trivial, whereas the p integration, using 


Eq. (B.2.12), is essentially 

o'6^~^\p')5^^^\x — Xi,) = {—T)^ ^ ^\x — xs)- (B.3.8) 


dp' 


Therefore the solution for (p, q) in the surface case can be simplified to 

where now 


xs = x^,{p' = 0,q' = 0) = 


1 + 1 /^ - gjp-^jy)) 

2y 


(B.3.9) 


(B.3.10) 
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B.4 The functions lE 


m 


In Eq. (B.2.12) we defined functions lE'^{xstar) given by one-dimensional integrals. These 


appear several times throughout this thesis. In this section we give explicit expressions for 
them. To simplify the expressions below, we shall commit the pre-factor 


VLd- 

(27r)^ 


4 \-l 

D-2 '^m 


(B.4.1) 


B.4.0.1 Even D 


If D is even and H — 4 > 2n, we can integrate by parts (H — 4)/2 — n times and get 

lE'^ix*) oc (-l)^-^a(^-"+™) [(l-x2)V+”^] 0(|x*|<l) (B.4.2) 


J X=Xi, 


= < 1 ), 


(B.4.3) 


where 




gi^-n+m) 


(l-x2)V+™ 


-H-m—n 


(i-x^r-i 




(B.4.4) 


k=0 


If Zl — 4 < 2n, which happens for {D = 4, n = 1, 2} and {D = 6, n = 2}, we use Eq. (B.1.7) 
but must proceed case by case. We get 


7g’^(x*) oc 7r6l(x* < —1) -|--arcsin(x*) 0(|x*| < I), 


12 


/f’^x*) oc -y^l - x^^dx*! < 1), 

/ 2 ’^(x*) oc 3x*i/l - xj^dx*! < 1), 

/o’^(x*) oc -x*7r6'(x* < -1) -h X* ( (arcsin(x*) - ^ ) -b y^l - xj ^dx*] < 1), 


TT 


I{ (x*) OC -^^{x* < -1) - 


^ ^ ^arcsin(x*) -b x*v^l - x^^ 


E’^ix*) ot (1 - x*)56'(|x*| < I), 

/o’^(x*) oc \ (arcsin(x*) -b x*a/ 1 - x^^ 

/^(x*) oc -(1 - x*)i6'(|x*| < 1), 

72’^(x*) oc 5x*( 1 - x*)i6»(|x*| < 1). 


k*| < 1), 


TT 


|a;*| < 1) -b < “1)’ 
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B.4.0.2 Odd D 


If D is odd and H — 5 > 2n > 0, we can also integrate by parts (H — 5)/2 — n times and get 
oc (-l)V-’^y’^ [(l-x2)V+”^] 5(|)(x-x*) (B.4.5) 


= (-1)- 


P^’”(x)5(2)(x-x*), 


(B.4.6) 


where 


Then we use 


[(1 - x^) V+™J = a 

k=0 

«(«(x) = Tj(-I)(x) = - - Ihiy 

dx 2yyr V xi ^ J 


(B.4.7) 


(B.4.8) 


The effect of the ^-function is to remove infinities (poles), so we can ignore it provided we 
retain only the finite part of the integral with the 0-function. So basically we have integrals 


of the form (finite part of) 


'-1 (x — X*) 2 


0(x — x^)dx. 


(B.4.9) 


We change variables to z = x — x* and expand the binomial, to get 


oc 


^+m-n k 


E «^E 


(B.4.10) 


X (1-x*)-^ 2 0(|x*| < 1)( (1 - x*)-’ 2-(-l-x*)-?' 2)0(x*<-l) 


If 0 < D—h < 2n, which happens for {D = 5, n = 1, 2}, {D = 6, n = 1, 2} and {D = 7,n = 2}, 
we define the negative-order derivative of (5-function as before and get 


oc 


D-5+m k k-j 

E <-^E( )7T^^ 


X \i-x^y+^-^e{\x^\<i) 

+ ((1 - x*)^+"-^ - (-1 - x*)^'+"-'^) 0(x* < -1)] , (B.4.11) 


where now the coefficients are defined by 


D—5+m 


5M(1_x2)^+«^ = ^ 


akx^. 


(B.4.12) 
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For some derivatives of the metric functions, some n = — 1 terms are actually required. For 


that case one obtains the same result as in Eq. (B.4.6), plus an extra term 


^(x) extra OC \P^’° - x) - P^’^ {-1 - x) 


m 
D-3 


(B.4.13) 


2^y TT 


9{l — x) 

(1 — a:)i 


+ (-!)- 


3 

( — 1 — X)2 


B.4.0.3 Derivatives 


We will also need to compute derivatives of the metric functions which, in the two-dimensional 
form, can be dg or dp (any p-derivative is trivial since the p-dependence has been factorised). 


Acting on the general expression of Eq. (B.3.3) or Eq. (B.3.9), there might be a term from the 


integration limits, and then there is the derivative of the integrand. Since all the dependence 
on p and q is through x*, and ^Im^^ = —Irn’'^~^, 

(B.4.14) 




In this way, it is straightforward to build the sources, namely for second-order calculation of 
Chapter 












Appendix C 


The asymptotic waveform 


This appendix contains details on the simplification of the integral formula for the metric 
functions in the limit r —>■ oo, including the factorisation of the angular dependence. It 
mainly supplements Chapter 

C.l The asymptotic limit in Fourier space 


Taking Eq. (B.2.10), we define the asymptotic waveform F(t,9), 


F{t,6) = lim rp 2 —F{u,v,p). 

r—^00 ar 


(C.1.1) 


and X as the asymptotic form of x*, 

1 


lim X* = — 


p' sin 9 


, 1 + cos 9 ,1 — cos 9 

-T + U --- \-V --- 


= z. 


(C.1.2) 


Noting that we get 


F{t,9) = / du' / dx' / dp' S{u ,v\p')I^’ 


dx 

dr 


(C.1.3) 


Then we dehne the Fourier transform of F{t, 9) with respect to the retarded time r, 


F{uj,9) = / dTF{T,9)e 


(C.1.4) 
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This integral is best done by inverting the function x(r) to t{x) and switching the integration 
variable, 


f+oo p-oo 


dr 




• +O0 


The X integration in Eq. (C.1.51 is essentially 


dx l!^’ s = ujp's\n6. 


(C.1.5) 


(C.1.6) 


From the dehnition in Eq. (B.2.12), this is 


dyAm(y)(l - l 8 ^°^^){y-x) 

(2tt) 2 j J-1 


O.D-4 

O.D-4 


D-2 

2 

-4 

D-2 




/ ^ d-2 

( 27 r) 2 


'-1 


"i 


{2tt)^ J-l 


where 


Ao = l, Xi = -{D-3)-\ X 2 = {D-3)-^{D-1) 


-It 


\-i 


(C.1.7) 

(C.1.8) 

(C.1.9) 

(C.1.10) 

(C.1.11) 


One recognizes in Eq. (C.1.10) the Poisson integral representation of the Bessel function of 
the first kind, 


Jv{x) — 


1 


dz{l- . 


y/TT [v+\) '^2/ 
Thus, inserting the volume of the n-sphere, 


(C.1.12) 


On — 


n + 1 

27r 2 

iW) 


we get 


Finally, 


dxl^’ ^{x)P^^ = i 2 +™sJd -4 . . 


(C.1.13) 


(C.1.14) 


C D — 2 poo D — 2 f 

F{uj,9) = —z“2 I dpp 2 jD^^^{ujpsm9)S ( oj 


1 + cos 9 1 — cos 9 


10 




, w- 


(C.1.15) 
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where 

S{x,y;p) = - J du J dv S{u,v, p). (C.1.16) 

S' is a double Fourier transform of the source with respect to u and v. The p integral can 
actually be cast as a (-D — 2)-dimensional Hankel transform, with argument casing. This is 
not surprising since the Hankel transform arises whenever a Fourier transform is made of a 
function with spherical symmetry. 


C.2 Factorisation of the angular dependence 


Now we shall see how the angular dependence factorises. Continuing from Eq. (C.1.15), we 
make the following transformations 


„ 1 . „ „ (s\.n9\^ "^l-cos^ 

p ^ R=-ujpsmO , to ^ O = to I I - - - 

^ \ A / Li 


and define a new transformed waveform 


such that 

We get 

f{0,Q) = - 


fiVL, 9) = ^/uJ'{n)F{uJ{n),9), 


du}\F{u},9)\‘^ = / dn\F{n, 


/ n x2 


^D - 3 Vsin6» 


1 — cos 0 \ 2 


X 0 


-1 


X / dRR 2 Jo-4 I ^(2R)S ( g; ^ 


1 + cos 9 1 — cos 9 2R 

.oj- 


(C.2.1) 


(C.2.2) 


(C.2.3) 


(C.2.4) 


(C.2.5) 


2 wsinyy^^^(f^) 

Now the CL symmetry comes in handy. From Chapter]^ we know that if F is 0(A:)J^ 

5(n,u,p) = (C.2.6) 


Therefore, and using dudv = p^dpdq, 


1 


S{x,y;p) = I dp I dq s^^\p,q) ^ 


—ixnn^ 2 {D 4 ) 


(C.2.7) 


Note that here F is either E or H so Nu = W = 0 in Eq. (5.1.161. 















124 


APPENDIX C. THE ASYMPTOTIC WAVEFORM 


Now we have 


p-2kiD-3) 

= (w sin, 

(C.2.8) 

y^{p) 

= Q^{R), 

(C.2.9) 

xp^-^ 

= n-^R^-^, 

(C.2.10) 



(C.2.11) 


Thus, 


, 1 + cos 9 1 — cos 9 

S ( u! --- ,(jj ---; 


2R 

uj sin 9 




2 {2R)^VD-3) 


and finally, 


f{n,9) = 


X 

X 


X ^ ^ X X 

y/D — 3 V1 — cos 9 ) V1 — cos 9 ) 

^ dRR^-‘^’^^^-^U^^j2R) X 

J dp J dqs^^\p,q)e~"'^^ in($(R)+^) _ 


(C.2.12) 


(C.2.13) 


At this point, the p, q coordinates are not necessarily advantageous since the R integration 
cannot be done as it is. So we might as well go back to 


u = qR 


D-2 


= pR-^^-^'> + ^R) , 


and rewrite Eq. (C.2.13) as 


-+m 


k-l 


-rtn a\ - * ^ V /l + cos^ 

’ “ /D^Vl-cos(9j Vl-cos^^ ^ 


X / dRR^jp-i j2R)s(n-^,n,R) 

Jo 2 


(C.2.14) 


(C.2.15) 

(C.2.16) 


This is Eq. (5.2.11). 




















Appendix D 


Gauge fixing of initial data 

This appendix contains details on the transformation of the initial data to de Bonder gauge. 
It mainly supplements Chapter]^ 


D.l General procedure 


Following the transformation in Eq. (3.4.5), we 


Eq. (3.4.8) is obeyed. The basic procedure is as 


need to find order by order, such that 
follows: 


1. Find the new de Bonder metric perturbations after the coordinate transformation in 
Eq. (3.4.5|) (£ denotes a Lie derivative) 






—)• 



—)• 



T Vfj.i' 

+ 2£^a)hW + (£ 


+ (a: - ^ G - 2) ^^ + • • • 


2 . 

3. 


Insert them in the gauge condition, Eq. (3.4.8). 


Use the wave equation to eliminate u-derivatives in favour of derivatives along the initial 
surface. 
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4. Now the gauge conditions are (schematically) in the form 


(D.1.2) 


where D is some differential operator and Cm ^ known function of (u, x*), which is com¬ 
puted at tt = O'*" from lower order gauge transformation vectors and derivatives 

of the initial conditions along the initial surface (no u-derivatives). 

5. Assume an ansatz 




(D.1.3) 


to solve (|D.1.2|) and find the functions and 




6. Finally, insert those in Eqs. (D.1.1), together with Eqs. (3.3.3), to obtain explicit 


expressions for the initial conditions in de Bonder gauge. 


D.2 First order 


The first-order metric perturbation transforms as 

/l(l) ^ /j(l) + £(1) +f(l) 

"'ai' ^ '"uu ' Su.iy I Si/,u • 


(D.2.1) 


Then the gauge condition in Eq. (3.4.8) becomes (dropping the superscript) 


— r.^,oiv 


(D.2.2) 


To eliminate the u-derivative, we use the wave equation. 


^uv — hya,ii 


(D.2.3) 


Finally, we get three equations for 


r)^,OiV 


(D.2.4) 


a = V yields an identity, a = j reads 


1 1 

‘^hyj^yy Lij ^y p ~hyj^-, 


0 . 


(D.2.5) 
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For a = u, we need to eliminate the u-derivative from the trace. Again from the wave 
equation, 


so we get 


h,uv — 4 ^ 1 ** ’ 


( ^u,v -if-) = i" ,22 H“ ‘^huu,vv hiu^iv “1“ 


(0-2KZ)-3) 


Therefore, the only component that changes is 

In particular, the radiative components do not change. 


(D.2.6) 


(D.2.7) 


With the ansatz in Eq. (D.1.3), one possible solution is = 0 and 

^(hi)=^(hi) = o, _CD). (d.2.8) 


(D.2.9) 


The term h{v, above makes it look like a second-order correction. However, when written 
in (p, g, p) coordinates, the power of p is the same as that of the original hul and agrees with 
Eq. ( 5.1.14| ) for A; = 1 and = 2,Ny = 0. 


D.3 Second order 


At second order, the first step is to compute the transformation induced by the previously 
computed From Eqs. (D.1.1), we must compute 


given the only non-vanishing components of 


(D.3.1) 


(D.3.2) 


The contribution from the first term is 




(D.3.3) 
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while that of the second is 




(D.3.4) 


All these vanish on it = 0, so there is no contribution from For the second-order gauge 
transformation, 

> (D.3.5) 

we proceed in a manner similar to the hrst-order case. The gauge condition is the same, 
Eq. (D.2.2), but Eq. (D.2.3) acquires a new term from the source 


Ihiice,uv — huct^ii + + T^v , 


(D.3.6) 


which also modifies Eq. (D.2.4) to 


- r)^,OV T ‘^^UOt,VV hia^iu p -|- „Tq, 


(D.3.7) 


Again, a = v yields an identity, hence we choose = 0 (and = 0 as before). Then 
a = j and a = u yield equations of the form 




,{ 2 , 1 ) 

'>j,VV 

= fuiv,X^), 


(D.3.8) 

(D.3.9) 


(2 1 ) (2 1 ) 

which can be integrated to find P ’ ^ and ^u’ '■ There is no need for the explicit solution, 


( 2 ) 

since the radiative components, /i) - , are unaffected. 





Appendix E 


Surface terms and late time tails 


This appendix contains details on the simplification of the computation of surface terms, i.e. 
the linear contribution from the initial data to the asymptotic metric functions, and on the 
late time behaviour of both surface and volume terms. It mainly supplements Chapter 


E.l Computation of surface terms 

In this section we specialise the results of Appendix to the linear contribution from the 
initial data. The source for the surface terms is 


S^^\u,v,p) = 25{u)d.F(-^\0,v,p) , F(^\0,v,p) = fi,^\p){v - <h(p))^0(u - cl>(p)), 

(E.1.1) 

where k = {1,2} and fl^\p) oc for E and H. Then, using the dehnition in 


Eq. (C.1.161, we get 


S{n-\n;R) = ■ 


(E.1.2) 


Inserting in Eq. (C.1.15), 


.fW(O) = fch 2^ (E.1.3) 

y U — 6 Jo 2 

As explained in Chapter we Eourier transform again to a new time coordinate t (see 
discussion in Section |6.2p, 

(E.l.4) 


1 


j:{t) = — / di}F{n)F 
2tx I 
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To perform this integral, we change the integration variable from R to T = ^ be- 

poo p^{oo) 


I dR^ I R = <I>-^{T). 

'o J^[o) 


(E.1.5) 


Then the fl integral is 




(E.1.6) 


yielding 


m = 


\/D — 3 y$(o) 


dT<^'{R)-^R'^Ji^,(2R)fl,'^\R)S^^-^\T - t) 


1 d 


,/D^ [<L'{R)dR 


k-i ( Ak) 


R —f^\R) 
^'{R) 




(E.1.7) 


where i? = $ ^{t). 


From the initial data in Eq. (3.3.3), we read 


EA): k = l, m = 2, fo{p) = - 


^'{p) 


EA): k = 2, m = 2, f^{p) =-{D - A) 


P 

^'{pf _ (11-4) 

p2{D-2) 


n -2 ’ 


4/9^ 


/fW: k = 2. m = 0, /„(p) = (D-3)^^ = tg54 


(E.1.8) 

(E.1.9) 

(E.1.10) 


and obtain 


Tui^mR)) 


— I 2 D-A 

— ■ R 2 JD (2R) , 

Vd^ -2^ 

D D — A D-6 
i 2 , =R 2 


/D^ 

i 2 y/D — 3 R 2 (^J D-4 (2R) + RJ D-2 (2R)^ 


(E.1.11) 

(E.1.12) 

(E.1.13) 


E.2 The contribution to the inelasticity 


The contribution of P{t) to the inelasticity is in Eq. (6.2.6). Recall that £‘^ oc (E P H)^. 


At = 1 there is only one term, eA)eA)_ The relevant integral is 

1 1 

i =D’ 


(E.2.1) 
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which converges for all D. 

At N = 2 there are two terms. In the relevant integral is 

1 2 1 
Jo ^ ziDT2 ■ 

This integral converges for D > 0. 

In the relevant integral is 

^ dR^JD{2R)(^Ji^{2R) + RJi^{2R)^ = 

This integral converges for D > 4. 

At N = 3 there are three terms. In the relevant integral is 

This integral converges for D > 0. 

In the relevant integral is 

I +i(2fl) {Jo^m + RJs^i2R)) = i (o _ p ^ 2) 

This integral converges for D > 4. 

Finally, in the relevant integral is 

Jo dRj^[jn^{2R)+RJD^{2R)'^ ^D{D-4){D-8) ' 
This integral converges for D > 8. 

These results are summarised in Table [621 


E.3 Late time tails 


In this section we compute the late time behaviour of the waveforms, both for 


volume terms. Following the strategy of Section C.l for the surface terms, we dc 
Fourier transform from fl to the new time coordinate t. The relevant integral is 


1 






1 


dLl e ^ . 


(E.2.2) 


(E.2.3) 


(E.2.4) 


(E.2.5) 


(E.2.6) 


surface and 
I the inverse 

(E.3.1) 
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where 


P = ^{R) +pR-^^-^^ - t, Q = qR 


.D-2 


(E.3.2) 


This integral can be made and equals 

i-‘^k-iQ{2k) (2\/X 

Therefore, we get 


X = PQ, 


(E.3.3) 


T{t,e)(x J dp J dqs^^\p,q)G'^{p,q;t), (E.3.4) 

where G^{p,q;t) is essentially the asymptotic (i.e. for an observer at null infinity) Green’s 


function in p, q coordinates and for the new time t, 

2k 


1 \ ZK pCSD ^ „ 

^ f dRR^- 




j^^^(2i2)0(X)Jo (2v/x) , (E.3.5) 


where 


X = PQ = qR 


D-2 


P + 


D-4 


R-iD-^)-A . 


(E.3.6) 


The action of the 2k derivatives distributes as 

2k 


j=o\3jj\X2 ^ ^ 

(“) {2VX 


(E.3.7) 
. (E.3.8) 


The R integral in (E.3.5) can be readily done for all but the very last term. So the contribution 
of the series inside the square brackets for G^{p, q; t) is 


2k-l 

E 


2k 


q-> 


1 d f j^^-2k{D-3)+j(D-2) 


J + \<I>'{T)dR 

where 

R = 


4* ^{t-p) , 


<L'{R) 


D = I 


Jd^+^{2R) 


t 


D-4 


-P 


(E.3.9) 


(E.3.10) 


p+ 1 , T» > 4 


Late times t >> 1 correspond to i? << 1. Thus we can study the late time tails by making 


a series expansion of Eq. (E.3.9) for i? ~ 0 
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The most relevant term is j = 0, since it has the highest powers of R (indeed, for k = 1 
this is the only surviving term because of the S{q) in s^^\p,q)). The waveform will only be 
square-integrable if its tail decays at least with Thus we conclude that E^^'^ and E^"^^ are 
square integrable for all D > 4, but is not: it grows exponentially in H = 4 and with 
t~°‘ for higher D. The condition a > 1 is only true for D > 8. 




Appendix F 


Asymptotic expansions for the 
Green’s function 


In this appendix we study the behaviour of the Green’s function for an observer at null 


infinity, when the source point is also near null infinity (Section F.l) and when it is close 


to the axis (Section F.2). This is important because a rapid growth or decay can induce 
numerical instabilities and loss of precision. It mainly supplements Chapter (the second- 
order calculation), but we present results valid to all orders. 


F.l Source point near null infinity 


When the source point is near null infinity, r/' —)• oo or, equivalently, Q' ^ {D — 2){D — 4) 
(for fixed Thus a convenient small parameter is 


,_ Q>-{D-2){D-A) 
^ 2{D - 2) 


(F.1.1) 


Then we have 


A'd-3£' = < 


I + 2e' - 2e' In Ae' + 2(^' - ^)e' + 0(e'2) , D = 4 
1 + {D - 4)1^^^ 


(F.1.2) 


. 1 + {D - 4)d-3^ 


+ 0{e') , 


D>4 
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and x* becomes 


X* = fo{y) - e'fi{y) + 0(e'^), 


(F.1.3) 


where 


fo{y) 


fiiy) 






l + f 
2y 


D = A 


1 + y2(^-3) 

2y^-3 ’ 

^-^' + ln4y 


y 


D>A 


-Ey^iD-3) + 

{D - A)yD-3 


D = A 


D>A 


(F.1.4) 


(F.1.5) 


and we have defined a short-hand 


„ _ 1 + {D- A)x>^i' 
1 + {D - 4)^^ 


(F.1.6) 


In the limit e' ^ 0 the integration domain collapses to the point y = 1 x* = 1. To show 
that the integral remains finite, let us define new variables w and rc*, 


y = 1 + Ye'w , X* = 1 - e'tc* 


(F.1.7) 


Then, expanding Eq. (F.1.3) around y = 1, 


w = 


m 


(/i(l) - w*) , 0 < u;* < /i(l) 


(F.1.8) 


This gives a finite domain for the integration (either in w or re*), and the shrinking measure 
dy cancels the divergent integrand. 


/ 



rsj 


y/pdWi, 

\/2uW 


0 ( 1 ), 


(F.1.9) 


Thus we conclude that there is no extra suppression factor in G^(A', Q') as r/' —)• cx). 
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F.2 Source point near the axis 


The axis p' = 0 corresponds to rj' = , so we take the limit —)■ r/' for fixed r]' < A 


convenient small parameter for this expansion is e = Q' 2 . Indeed, in D = A, 

c/ / o ^ 1 + Q' ~ 1 

' Q' ViW + i 

~ -4e + 0(e^), 


(F.2.1) 

(F.2.2) 


and, in H > 4, 


1 + {D - 
1 + {D - A)^p' 


r-D + 3 / Q' 

F + D-3 VQ' + 2(l-^ 


~ 1_ 2P-2)(D-4) 

“ (D - 3) £ + F 


(F.2.3) 


(F.2.4) 


Then we solve ^'{P'lQ') for P'{^',Q') and substitute in A', 


- 21n- 


- 2(C - r), 


T» = 4 


1 + {D-A)T^^' 
1 + {D - 


D-2 

q-{D-3) 


D-A 


, D>A 


(F.2.5) 


Expanding, we find, in Z1 = 4, 


A' = 1 - 2(^ - p) + ln(4e2) - 4e + 0 {£^), 


(F.2.6) 


and, in H > 4, 


where 




(F.2.7) 


1 + {D -A)T^^ (A*-3) 


(F.2.8) 


Then we scale the integration variable y. 


D = A 


D -A\ d-3 


z , D > A 


(F.2.9) 
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such that Eq. (7.5.6) becomes 


= fo{z) + efi{z) + Oie"^), 


(F.2.10) 


where 


Mz) = < 


—I — 2{^ — T]') + z"^ — 2\n z 

Tz ’ 

{D - 4)z2(^-3) - 2{D - + (^ - 2) 


2{D-2){D-A)z 


D-3 


Mz) = < 


D = A 


, D>4 

(z 


D = 4 


D>4 


. (F.2.11) 


(F.2.12) 


F.2.1 Inside the light cone 

Let us start with the case For a consistent, well-defined solution in the limit e —>■ 0, 

we define a new integration variable w through 


z = zq + ew , 


(F.2.13) 


(F.2.14) 


where zq solves fo{zo) = 0. Then the equation for x* becomes 

X* = fo{zo)w + fi{zo) + 0{e), 

which we can invert for t(;(x*). 

Then, in the integral which defines the asymptotic Green’s function, 

D — 4 

Gl(A',Q')oc / dyrC’\x.), a = - {D - 2,){2k + N^ - N,), (F.2.15) 


we substitute Eqs. (F.2.9), (F.2.13) and (F.2.14) and transform to an integral in x* 


G^{A', Q') J dxi, /^’°(x*) (1 + eu;(x*))“ + ... , (F.2.16) 


where 


/3= < 


— (a + 1), D — 4 
(y P I 


(F.2.17) 


D-3 


D>4 
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Now observe that (see Eq. (B.2.12)), 




(F.2.18) 


vanishes for n < (E> — 4)/2 + If each extra power of x* comes with another power of e, 
then the first non-zero term will have a power of 


^ E»-4 

/3 H---h m . 


(F.2.19) 


F.2.2 On the light cone 


When 7]' = ^, we have 


zo = l, /o(l) = /^(l) = 0, 


(F.2.20) 


so we require a different definition for w, similar to Eq. (F.1.7), 


z = zo + y/ew . 


(F.2.21) 


We get 


=-(x* -|- 1). 


Thus, in Eq. (F.2.15), the transformation from y to w yields a power of e equal to 




(F.2.22) 


(F.2.23) 


For the second-order calculation of Chapter this is stronger than outside the light cone. 


This is obvious in even D but by analytic continuation is also true for odd D. 
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